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A MUDLL  roK  THL  DOPPLER  SPREAD  GF  DACKSCATTERED  SOUND  FROM 
A COMPOSJTE  ROUGHNESS  SEA  SURFACE 

by 

Herwai’d  Scliwaiv.e 


ABSTRACT 


A tlieoret  ical  model  t'oi'  tlie  dopploi’  spread  of  bacRscal t ered 
acoustic  waves  fi'om  the  rough  sea  surface  is  dcsci'ibed  on  llie 
basis  of  a two  coniponent  st  I'ucturo  (facet  model)  of  t b.e  sea 
surface.  The  facet  statistics  are  derived  for  a general  gaussian 
sea  surface.  The  results  are  evaluatfsd  for  a Pierson~MosKowi  tr. 


spectrum.  A proceduie  for  choosing  the  facei  length  is  developed. 
Approximate,  simple  i'ormulae  foi  the  Fiorson-Moskowi t z case  are 
given.  The  general  result  for  the  doppler  sp’’ead  of  the  back- 
scattered  acoustic  wave  is  evaluated  for  a sea  sui'face  spectrum 


due  to 
f erent 

f I"*  rn  1 


Scot  t and 
wi ndspeeds 

! 4 C ^ O S n 


examples  foi'  dopplei'  spectra  are 
, wind  directi oijs,  grazing  angles 


given  for  dif- 
and  acoustic 


INTRODUCTION 


The  important  pioblem  ot'  scattei'infi  of  acovistic  waves  fi’om  the  lough 
nio\  i nj.’,  sea  surface  has  olten  be<’:r.  i nvest  i j;at ed,  but  a geuci'al.  exact, 
and  tractable  solution  has  nut  e<  been  oblained, 

lor  special  cases,  solutions  aie  1 .lowii,  e.f;.  foi’  the  cases  where  the 
amplitudes  of  the  sea  sui’i'ace  waves  are  much  smallei’  [Ref,  ll  oi’ 
much  lai’f^ei  [Ref.  2]  than  the  length  of  the  incident  wave, 
lor  arbivraiy  i oughnesses  of  the  sea  surface,  approximate  solutions 
a>  e ubtaLiieU  by  applying  the  small-scale  bacbscat  t er  i ng  results  to 
a composite-roughness  sea-surface  model  [Refs.  3,  4,  5,  t'J.  This 
means  that  the  small  waves  with  lengths  up  to  some  10  cm,  causing  the 
1 esonant  oi  "Bragg"  scattering,  are  cari ied  hy  long  waves.  The  long 
waves  ai  c approximated  by  plane  faeots  whose  movements  depend  on  the 
large--scale  roughness  of  the  sea  sui’facc. 

In  this  paiioi  , the  interest  is  focussed  on  the  frequency  spreading 
ol  a backscatt ci'cd  monochromatic  acoustic  plane  wave.  This  case  is 
ol'  special  interest  in  active  sonar  applications  because  the  moving 
sea  sui'fai:c  limits  the  detectability  ol  lilovily-moving  targets, 
opecial  effort  has  been  devoted  to  obtain  a set  of  formulae  that  are 
simple  to  use',  in  already  existing  sonai’  models,  especially  in  the 
RAIBAC  program  system  developed  at  SACL.ANTCEN  [Ref.  7j. 

The  main  text  of  this  paper  contains  the  lesnlts  obtained,  logethei’ 
with  some  physical  explanations.  All  mathematical  derivations  and 
proofs  are  given  in  appiuidiccs,  ' 

In  t he  III  lin  text,  the  (existing  resonant  scattering  theory  is  first 
icviowed  and  the  facets  are  then  introduced.  The  statist,!  cal  prop- 
erties c>f  the  facets  are  given  in  terms  of  the  sea-snrface  covariance 
function  or  the  t liree-dimensi  onal  sea-surface  spectrum.  The  concept 
is  applied  to  the  small-scale  roughness  I'esults,  Seteral  special 
cases  ai'e  considered. 

The  geiuu’al  formalism  is  worked  out  in  several  examples  of  idealized 
sea-surf  ace  sjiect  ra  and  the  computer  i'esults  arc  shown. 


1 . RESONANT  BACKSCATTERTNG  THEORY  FOR  SNiALl,  SEA-SURl  ACE  ROUGHNESS 


Tt  U'  I'v  KhnWM  in  F i _ 1 is  n.snd.  A S ( Xn^. . 

Illuminati's  a sui’face  ai’ea  A = ab  witli  a plane  wav'e  of  fi'equency  fo 
aiid  pressure  amplitude  p^,, , The  distance  Tq  is  considered  inucli 
larger  than  »/ A . Without  loss  of  generality  tlie  y— coordinate  of  S 
IS  set  to  zero.  Besides  the  specular  reflection,  the  lougb  sea  sur- 
face causes  a scattered  field  in  all  diiections.  F'^r  backscattei'ing, 
soui or  S and  receiver  R coincide. 

The  doppler  spectayun  $(f)  of  the  ha  ckscat  t erod  wave 
Fourioi'  transform  of  its  correlation  function,  viz: 

OL 

J Eip^(t)p(t+T) } expj-2n  jfx!  dt 

.00 

where  Ei,.,|  denotes  expectation  and  ^ the  conjugate  complex. 


p(t  ) is  the 

[Eq.  1] 


A iiurmal  . ,.»i  ion  iy  introduced,  yioJdiiif;  a backseat  1 or iny  dopplci' 
dt'ns  i t V 


cp(l')  = 


Mr)rg 

Po  A 


[Lq.  2] 


Startiny  Ci'oiu  Ref.  8,  for  example,  the  basic  oqnzi.ion  t'oi'  the  back- 
scat  tor  Liiy  dopplei’  flensity  t'oi’  the  slightly  rciiyb  sea  surtace  is 


cp(t')  = (4n)‘"  sin*  Yol<oXa  (-2kocos  Yo  > 0>  I'-iy  ) • 

1 fn 


[bq.  3] 


In  bq,  3,  ko  = - — \=  — is  tbo  wave  parameter  of  the  acoustic  plane 

Ao  C 

wave  and  ^a(k^,k  ,f)  is  t!ie  f hi’ee— dime.nsiona  1 wave  pai’aniet ei’ 

fi'equency  spect  I’uin  of  the  sea  surface.  The  wave  parameter  k = ^ 

was  inti’oduced  hei'e  instead  of  the  wavi^  numb(U’  3:11  , because  it 

\ 

offei's  some  notational  ad\'antage.  Equation  3 is  dei’i  ved  in  detail 
in  App.  D.  When  integrated  over  f,  Eq.  3 becomes  the  back- 
scat  tci’ing  strength,  known  from  the  literatui'e  [|Reis.  4]. 

The  general  small-scale  roughness  result  is  simplified  b>'  the 
assumpt  ion  that  the  dispei'sion  I’elat  ion 

- gk  [Eq.  4] 

is  valid,  where  g is  the  gi'avity  constant. 

Using  the  lelat ions  in  App.  A,  Eq.  3 is  written  in  terms  of  the 
frequency  angle  spectrum  Fs(f,  cp) 


2g"  k(,sin*Yo 


cp(f)  = 


1, 0) 


[bq-  5] 


[W(f/,n)  Mf-fo-l'o*) +W(bo',  0)  6(i'-fo  + bo^)], 


where 


1 o ~ V 


IT 


[bq.  6] 


is  the  dispersion  relation.  In  Eq , 5,  b?  ( i j w)  snd  W(f,  cp) 
ar<'  the  frequenc>'  angle  spectrum  and  the  wave-mixing  fnnet  iort, 
iW'S[)ect  i vely  . 

The  mixing  funct  ion  describes  the  I’at  io  of  incoming  and  outgoing 
sui'face  waves  for  a given  wave  paramet  ei’  vector  oi’  a given  frequency 
and  direction.  The  introduction  of  tliis  function  leads  to  the  same 
lesults  as  the  theory  of  Neumann  and  Pierson  [Ref.  10,  p.  330]  but 
it  is  felt  that  the  formulation  used  heie  is  sinrjler.  I For  details 
of  the  definition  of  W and  its  iclation  to  the  wave  paramet  tu’ 
frequency  spectrum  X3  the  I’eader  is  I'cferi’cd  to  App.  A. 


4 


A sinipTr  pi;»ctical  catiij  is  the  omii ' <1 1 reel  i ona  1 sea— stu'face  spccli’uni 
VI  111  W = 0,5,  which  seenis  i t^asmiah  n tor  short  sui’face  wa\es.  In 
tills  case  oiK^  ohtaitis 

-I'o-to  ) t-  6(t-to  + to)]. 

[liq.  7] 


cp(t)  = 


siii*Yo^4 


2tti 


.v^ 


l‘l  (*0 


Lqnatiems  5 ami  7 are  coiisidei’od  to  be  valid  midei’  tlie  following 
Coiidi  t tons : 


1) 


2) 


3) 


H I 


Xq  *^h  * 2TTsiii  Yg 

ajj!  si  amlard  del  iat  ion 
of  the  sea  sni’facc 
ht'ighl 

X «7a 

0 


,^»/A 

— >0,05  III 

Iv 


"Bragg"  scattering, 
RavJejgh  theory 
[Kef.  51] 


implicit  assumption 
I'or  I’qs . 5 and  7, 

e'xplanat  i on 
set!  App.  B 

Assumption  of  incidemt 
jilane  wave 

Approxi  mat  i tin  I'li  1 e for 
tile  validity  of  llic 
dispt.'i’sion  I’e'l  at  ion, 


[tq.  8] 


[Kq.  9] 


[Eq.  10] 
[Eq.  11] 


As  most  acoustic  wai  t'l  engt  tis  in  sonar  applications  art>  be'twoen  0,1 
and  1 111,  t ht'  condition  of  Eq,  8 is  of  t tMi  not  fullilled.  Therefore 
it  j ,s  nect'ssai  > to  introduce'  tin'  composite'  roughnt'ss  st!a  surface! 
mode'l,  wliieh  c ii  cunnt'iit  s t h es  ctnitli  t ion , Th  i .s  modt'l  is  devt'loped 
in  tht'  next  chapteir. 


2 . nOPPIER  SPECTRUM  FOK  A COMPOSITE-ROUGHNESS  SEA  5URKACE 

2 , I I nt  roduct  ory 

T1  le  dciivations  of  this  chapter  ai'e'  based  on  t lie  facet  model  as  used 
by  Bachmann  [Rt'fs.  5,  14]  foi’  tht'  calculation  ot  the  backseat  t ei' j ng 

st  rt'i.gt  h , 

The'  facet  mode]  assumes  basically  a t wo-eompone'nt  structure!  of  t lit! 
rough— mo  \' i ng  se'a  surface:  the  siual  1 — seal  »'  lu  gh-f  requenc>'  ’'ouglmess 

(.lippleis),  rt.'sponsi  bl  e 1 oi'  the'  I’e'seinant  scatte.iing  are  cat'  'ed  by 
tht'  1 e)w— f re'qut'iic>-  large'— scale  se'a— surf  act'  waves  (swell).  Those' 
cari  iei  waves  art'  locally  appi'oximat  I'd  by  jilane  facets  of  finite' 
e'Xtt!nd.  The  HKivement  s of  these!  face'ts  ha\b  influence's  on  the'  back— 
seiatee'iing  slreuigth  and  the  tlopplei'  shifts  of  tht!  backseat tei'od  sound. 

In  the  first  part  of  1 li  i s chapter  the  1 act't  mode'l  is  tlcee' 1 tiped  and 
the  st  at  ist  it!S  ol  the  face'ts  are!  calculate!d.  In  the  se'cond  part,  the! 
t'ltt'c'ts  of  the  t ace't  moNe'iiu'iit  s ein  the-  small-scale!  roughnt'ss  re'sults  ai'c 
C'ons  1 dert'd , 


5 


2.2 


1 acc(  staliwt  ii's 


The  I’acet  ai  e st  ai  is!  ica  1 ly  described  by  a five-dimensional  random 
j>i’ocess:  i Ik;  inclinations  Cx  and  Cy  and  the  ^’elocities  Ux*i  Uy, 

and  U.T  . If  the  iticidc'nt  a const  ic— wave  vector  lies  in  flu;  x— /;  plane, 
the  influences  of  Gy  and  u^,  can  be  considei’ed  small  and  are 
neglected  [Ref.  14].'  Consequently  one  deals  with  a three-dimensional 
procK'ss  onl\  and  the  facet  can  be  I’opresent  ed  by  a sti'aight  line, 
thus  reducing  the  mat  h.emat  i ca  1 effort  . 

As  it  is  assumed  that  the  sea  sni'fac<*  is  described  as  a gaussian 
pi  ocess,  the  facet  ino\enu;nt  is  also  coiisid«;red  to  be  normal  disti’ib- 
uted.  As  i 1k'  mean  \alo«.-s  of  tlie  stochastic  variables  e = Gx(f)j 
ux(l)  and  u^  ( t ) are  zero  by  definition,  tbeii’  covariance  matiix 
is  sufficient  for  theii’  complete  description. 

The  gcoiu'ti'y  sliown  in  I'ig.  2 is  used  in  the  calculation.  To  obtain 

an  expression  for  C(t)  and  b(t  ),  t lu*  sea  surface  tunction 

h(x,  y,  t ) I ^,=0  is  approximated  by  a 1 <-ast  — sqnai'i;-er  ror  straight  line 

f(x,t)=c(t)x  + b [Eq.l2] 


for  t lie  horizontal  I'elocity  Ux  t I'e  following  argument  applies. 

Ti’om  the  linear  t heoiys'  ol'  sui’face  waves  it  follows  that  the 
horizontal  particle  lU'locity  u in  t !ic  x-diiection  is 

iii  - -g  y,  t ) . [Eq.  16] 

9t  9x 


b 


7 


is  defined  as  th<'  mean 


The  lioi  J ontal  velocity  of  the  facet 
hoiizontal  particle  voiocity  on  the  facet 


[Eq.  17 


Insei’t. ing  Eq.  3 6 into  Eq.  3 7 gives 
t 

~ ~ |^[h(|-  , y,  l.^)-h(-'|.,  y,  ti)J  dt^  + u^(to  ) . [Eq.  18] 

t 

0 


u (t) 


I. 

7 

1 r H 

— , udx 


L J. 


From  Eqs.  13,  3 5 and  17  the  variances  and  ci’oss- variances  of  the 
facet  movement  are  calculated.  The  results  are  given  below,  the 
derivation  of  the  foiniulae  is  found  in  App.  C 


and 


E*  u‘-  ! 
t.  X 


L 

M r 

1= 

o 

L 


C(x,  y,  t) 


3x  ^ 2x3 


+ — — I dx 


2 ra=c(x,  y,  T)  ('i.i)dx, 

L ,1^  5t  l; 


[Eq.  lb] 


[Eq.  20] 


T -U 


,,2 

[f  [2C(0,  y,  t)  - C(L,  y,  t)  - C(-L,  y,  t)]dt  dti  . 

^ y=0 

[Eq.  21] 


The  cross  variances  arc 


Eje  u^ 


and 

{ e 1 


12  P ^C(x,  y,  T ) 


J 


9t 


o 


— - 

|y=o'  1 '' 


dx 


T=0 


E { u u I = 0 . 
t.  ' X z ’ 


[Eq.  22] 


[Eq.  23] 


Using  the  results  of  App,  A,  the  above  equations  can  be  written  in 
terras  of  the  wave  parameter  frequency  spectrum  or  the  frequency  angle 
spectrum.  As  the  latter  is  of  importance  in  the  applications,  the 
formulation  of  the  facet  statistics  in  terms  of  p2(f,  cp)  is  given 
here , 


S 


It  is 


+ TT  «• 


i e' 


144  r r ,,  ai.  n • aI  -»3  ,r^.. 

— — t I — — (La  cos  "^-2sin-^J  ciLdcp, 


J J 

-no 


[Eq.  24] 


+ n 


K la"  I - — 

Z ' ,3 


2g  p p lg(f,  fQ)(l  — ct)saL)  ...  , 

J J ; ; didCp, 


-no 


( 2nl'  cos  cp)‘ 


[Eq.  25] 


E,  ju"  ] - 

t X ' 


+ n “ 

p p FgCf,  cp)  (1  - cos  aL) 
~7T  I i 

L U L- 
- TT  O 


(2nt-)" 


dt'dcp  , 


[Eq.  26] 


and 


E^  1 s 1 


i-n  « 

— i r I-’s  (t,  cp)  (2W(f,  cp)  - 1)  2nf 

j 3 J i.J 

-n  a 


sin  aL  2 

+ ( 1 - cos  aL  ) 


.3 


dfdcp  , 


ca 


[Eq.  27] 


where 


(2nf)"  cos  cp 


[ Eq , 28) 


The  dei'ivation  of  the  above  I’oi'mulac  is  also  found  in  App.  C, 

The  horizontal  v«'locity  is  uncorrelated  to  C and  , as  is 

seen  froni  Eq.  23.  As  C , u.^  and  ii^  ai’c  assumed  to  be  normally 
dist,ri  but  edj  u^  is  thei'cfoi  e independent  of  e and  u^ , Then 
the  tliree— dimens ional  probability  density  Wg(e,  u^,  u^^.)  has  the 
form 


W3  (e,  Oj,,  n^)  = wj  (c. 


wi(u^) 


where 


w,  ( u ) ^ 

i.  \ V / 

J 2n  CT 


exp 


N(u^,  0,  a^) 


[Eq.  29] 


[Eq.  30] 
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exp 


^3(2,  'ij.) 


] 


2TT(T^a 


z;  € 


J) 


[Eq.  31] 


Tlic  abbre\  lal  ious  in  Eqs.  3b  and  31  ai'c 

^ ^7.^  v^Et  f,  =-  y E^  1 [Eq.  32] 

ai'd 

E i C u I 

P =-=  ^ . [Eq.  3 3] 

cj„  a,„ 


Equation  2S^  doscribet;  i ho  I’accl  .statistics  and  i.s  used  in  the 
t'<;Llo\  ing  lor-  the  calculation  ot  tlu^  dopploi'  density. 


2 , 3 Choice  of  the  facet  length 

Tlie  fac;et  statistics  as  derived  in  the  preceding  section  depend 
entirely  on  the?  racet  length  L.  In  practical  cases,  this  L has 
to  be  chosen.  The  choice  is  dii'ecled  b\'  two  ini'!  nonces , Firstly, 
the  facet  length  should  be  great ei'  t.hatt  the  wave  length  of  the 
incident  acoustic  wave  to  reduce  the  finite  aperture  error  as  much 
as  possil)le.  That  lt^ads  tcj  ilit  condition 


L » Xo  , 

[Eq.  34] 

Secondly,  t ]io  sea-stu'face  roughness  , -71^(1)  on 
much  smaller  than  t lui  acoustic  wave  leiigtli 
Rayleigli  condition  [Eq.  8j 

tlie  facet,  should  be 
to  fulfil  the 

Xo»aj.(L)  . 

[Eq.  35] 

It  seems  to  be  cxtrcmel}.-  complicated  to  calculate  the  finite 
apci-ture  error  and  tlie  finite  roughne.ss  erioi'  and  to  choose  that 
facet  lengt  h,  Lopt  ;>  where  the  sum  of  bot  h is  a minimum.  Therefore 
the  following  procedure  for  choosing  the  facet  length  is  performed: 

A ruimher  N is  cliosen  that,  transforms  the  two  aboie  inequalities 
[Eq.s.  34  and  3 5]  to  the  t^quations 

I,  - Xo  " N 

[Eq.  36] 

an  d 

[Eq.  37] 

Combining  them  gives 


L • a.  (L)  = 


[Eq.  38] 


Tills  implicit  equation  has  to  be  solved  to  find  the  facet  length  L. 
The  number  N is  a quality  nuinbei’  tiiat  shows  how  far  L is  away  from 
the  limit  whei’e  the  conditions  for  tlic  validity  of  the  composite  sea 
surface  model  are  no  longer  valid.  If  N is  less  than  about  1,  tliis 
indicates  that  the  facet  model  will,  yield  no  correct  results. 

To  solve  Eq.  38^  an  expression  for  a.  (l)  has  to  be  found.  The 
roughness  r(t)  or.  the  facet  is  ^ 


r(x,  t)-  h(x,  y,  t)l^^.Q-f(x,  t)  , 


[t:q.  39] 


wlierc  f(x,  t)  is  the  facet  equation  in  Eq,  12. 

TliC  mean  .sqnai  e ei'roi'  p*"  (t  ) is  a I'andom  variable  in  time 


] r 2 . . . 

L J r (x,  t;  dx 


[Eq.  40] 


The  sea— surface  roughness 


Cj.^(t)  is  the  expectation  of  p"(t) 


E„  1 P^  { 


[Eq.  41] 


tsiiig  Apf).  ky  is  it  shown  in  App.  C that  this  takes  the  form 

1 

E^ip"  i = C(o,  o,  o)  j^C(x,  y,  t)  --p  + ^)dx  [Eq.  42] 

o 

or,  written  in  terms  of  the  frequency  angle  spectrum  p2(f,  cp)  , 


+1T  CO 

E!pMt)}  = JJ  F, 


COS  at — ^ sin  aL+2H — 

at  1“  a 


df  dep  , 


[Eq.  43] 


where  the  factor  a is  seen  from  Eq.  28, 

With  Eq.  42  or  Eq.  43  it  is  possible  to  solve  Eq,  38  by  an  interactive 
procedure  for  a practical  case,  where  C(x,  y,  x)  or  Es(fj  tp)  is 
loiown  explicitly.  This  will  be  inA'^est  igated  in  Chapter  3. 
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2.4  The  Doppler  spectrum 

The  doppler  spectrum  for  an  arbitrarily  rough  surface  takes  into 
account  the  facet  movements  whose  statistics  were  described  in  the 
preceding  chapter. 

The  gra^.ing  angle  Yo  incident  frequency  f^  are  modulated 

by  the  facet  inclination  e(t)  and  the  facet  velocities  Ux(t)  and 
*iz;(i)  respectively  and  accordingly  the  small  scale  result  in  Eq.  5 
is  changed. 

The  modulation  changes  Yo  ^ind  t'o  to  statistic  variables  of  the 
following  form 


1)  Influence  of  the  facet  slope 


Yf (t)  = Yo  + e(t) 


2)  Infliience  of  the  facet  velocities 


fp(t)  f^  ■+ [u^(t  )sin  Y£-(t)  - u^(t)  cos  Yp(t)] 


To  obtain  the  backscatt  ered  doppler  density  (p(f),  the  weighted 
sum  for  all  C , Uj^,  and  u^  has  to  be  taken.  Using  Eq.  6 one 
f i nds 


pen 


W(f^:  , rr)  5(f  - - f^)  + W(f^  , o)  5 ( f - f £.  + f^,)  ] , 


w,  ( c , 11  , u ) dedu  du 

^ ’ z ’ s.'  z X 


This  integral  can  be  solved  in  closed  form  for  a small  facet 
slope  c and  other  minor  neglect  ions  that  ai’e  found,  together 
with  a derivation  of  the  general  result,  in  App.  D.  The  general 
formula  for  the  dopplei’  density  tp(f)  is 


cp(f) 


4fo  sin  Yo 


[/(f-fo)®  Fs  (f*  , o)  W(f*,  tt)  I_ 


+ y(f+f‘o)®  Fj.(f  , , o)  W(f.  , o)  I,]  , 
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f 


wher  V, 


H cos  Yo  ^'o  ) ® 


^ i *^1  ) i 3 (1  - Pg^)‘ 


■ + <,  + ) 


3a|  + 6ag  m^  + m*  ^]  . 


In  Lq,  4P  the  abbreviations  arc 


(f  - fo  + t'o  ) c 
2 to  s -in  Yo 


J C^.  + C^.  t an^  ' 


tan  Yo 


® C'^x  Pc  X 


+ CT^.  tan^  Yo  Oi  tan  Yo 


CT,.  P, 


u_^^tan^Yo 


ul  + a-;  tan*'  Yo 


— + tan  Yo  • 


From  this  general  result  evei’y  I’equired  simplification  can  be 
easily  deduced.  For  example,  if  the  facet  slojie  is  considered  of 
little  influence  in  a cei’tain  case,  the  expi'ession  1 in  Eq.  4P 
reduces  to 


T_  - lv(u_  , o , ai  ) tan*Yo  • 


[Eq.  54] 


Ano\,lHjr  special  ease  is  the  assuniptioii  of  an  onuiiUiecci ional  sea 
sui'face  spccti'iim  and  a mixing  function  W = 0.5  (seeApp.  A). 

Using  Eq.  A.  26,  the  dopplcr  density  foi’  this  case  can  easily  be  found. 

The  geiicra  l I’csult  of  Eq.  47  has  l>ceii  programmed  on  a compiit  er  <'or 
sea— surface  spec! '-a  of  pi'actical  importance.  These  applications  arc 
disenssed  in  the  ne.:t  chaptci’. 


3.  APPLICATIONS 

3.1  Examples  of  sea— surface  spccti'a 

Foi'  demonstrating  the  use  of  tlie  mode]  described  in  the  pi’eceding 
chapi  ers  a sea— surface  spect  I'um  needs  to  be  chosen  as  an  example. 

This  sped  I’lim  should  have  a well  defined  frequency  and  directionality 
variation  with  the  wind  speed  and  wind  dii-ection.  On  the  other  hand 
it  should  be  mathematically  simple  enough  to  allow  the  calculation  of 
the  facet  st  atistics  with  not  t.oo  much  ef i'ort . Therefore  t-he  well 
known  Piei'son— Moskowit  z specti’um  was  chosen  foi’  t lie  facet  statistics 
in  slightly  modified  foi'm.  It  can  be  found  in  App.  A,  Eqs . A,  50  and 
A.51>  and  is  igiven  here-  for  convenience 


Fa  (F,  «P)  Fi  (f) 


COs'\cp  - ePo  ) 


n 


where 


Fi  (f) 


ag~ 


(2tt)M’5 

0 , else 


^ 4/ 


and 


TT 

cos*^  cp  dep  , a = 0.0081  and  g = 0,74. 
-n 


[Eq.  55] 


[Eq.  56] 


The  mixing  function  W used  for  these  calculations  is  given  by 
Eq,  A. 40. 


For  the  calculation 
w'as  employed,  which 
of  the  sui’face— wave 
is  found  in  Eqs.  60 


of  the  doppler  density  a more  flexible  spectrum 
takes  into  consideration  the  fi’equeucy  dependence 
dii’ecti vity , This  spectrum,  due  to  Scot^  [Ref,  IS] 
and  A, 64  of  Appendix  A. 


3 . 2 Variances  of  facet  slope  and  velocities 

In  this  section  the  geneial  formulae  for  the  facet,  statistics 
[Eqs.  24  to  28]  are  evaluat  ed  for  tlie  Pierson-Moskowitz  sea-surface 
spectrum.  Again  the  derivations  are  found  in  App.  C,  the  results  are 
given  in  the  following. 
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The  staiHi„e  |.oinl  is  tl.c  P i,-.  so..->'o»ko..iis’ 

Vn  This  snoctrulu  is  insoi'tod  xiito  bqs . -4 

ov?r  ilK-  fr.M«o..cy  f is  ,.error,».,d,  yielding 


a V 


J 


p „ \ 18  48 

•'  cos”(®  - %)  ^ 


Wq 


u; 


U, 


+ cos  Uq 

\ 

- Ci( |uo  I ) I ) 


n 


2a^’“ 


CT„  - 


r cos'^  ( tp  - CPo  ) 1 1 ^0 


V 8 ■ o 


1 


- V) 

Uo  / 


sin  ug  Iuq  } I \ I 

+ — ^ Si(  1 Uq  I ) > dcp  , 


X 


6 0 ( ^0  \ 
J cos"  ( Cp  - tpo  ) I 4 + cos  Uo  \^-l+  — j 


\ 


Uo  sin  Uo  )Uo  I t ^ i , 

+ — Si(  Uo})V  dcp  , 

2 2 ' 


and 


E,  1 c u 1 = o 

‘ /.  t /:  Q 0* 


M + — . 

. '''  2 

12(]-2Wo)av*  p „ 

■'  cos  iCp  - ’ * 


! 


J 


L CP  -H 


si  n Ug  I — g-  + l<3  ) i‘  cos  Uq 


4^3 


Ut,  Uo 


("I 

-21-5  TfTO 


y dtp , 


ks 


_ + _ + k5  Uo  ) - ^^3 


where 

Uo  = ao  cos  cp 

and 


(ZTTfy)^  L ^/T  Lg 


g 


the  iorni  of 
Then  the 

[hq.  57] 

|]Eq.  S8] 
[Eq.  59] 

[Eq.  60] 

[Eq.  6l] 
[Eq.  62] 
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The  constanl  s kj^  to  ks  aro  k^  = - 10/63^  ha  = 8/l89> 
kj  — - k^  = — 4/63j  and  kj  = IO/I89.  The  luiictioiis  Si,  Ci 

and  S af<3  The  sine  iiiicfejeal,  the  cosine  iiite^jral  and  the  Fresnel 
intt^gral  [Rei,  13]>  respectively. 

Tlie  ^•al'iances  in  Fqs,  57  in  hO  liavc  lu;eii  prograiiuucd  on  a digital 
coiuputei’,  using  nuiiiei'ical  i nl  egi’al  i on  techniques.  lnsp<;clion  of 
flu  equat  ions  shows  that  the  va'’iances  ai’o  I'uuctions  of  llu; 
f o 1 1 ow  i ti g f o 1‘ni 


Zib  f Thq.  63] 


V 

ai\d 


Thei'efore,  the  results  ai'e  shown  as  a funct  ioti  of  L/v^  rather 
than  the  facet  length  L,  thus  eliminating  the  variat^ion  of  the 
results  as  a funct  ion  of  the  windspeed  . 

The  remaining  pai'amot  ers  ai’O  the  direction  angle  CPq  , tlie 
direct  ion  angle  cp  of  t he  mixing  funivt  ion,  and  the  exponent  n 
of  t l\o  cosine  dii’cct  ivity  funct.  ion.  In  Figs.  3 t o 12  t he  facet 
statistics  are  shown  tor  epo  = = 0®,  45®  and  90®  and  n — 2 

and  n ~ 4.  Fi’om  the  figures,  the  facet  statistics  can  be  found 
for  pi'aotical  cases, 

Tlur  impoi’tant.  I'egion  for  L/ v®  between  0,1  and  1 s^/m  is  shown  in 
a linear  scale  in  addit  ion  to  the  genoi  al  doxible  logarithmic  plots. 
Ihe  facet  slope  divergt:s  for  L Oj  therefore  the  curves  in 

figures  3b  and  7h  start  at  the  maximum  value  on  tlie  a —scale. 

For  < 0.1,  Figs,  3a  and  7a  should  bo  used, 
v'' 


F^or  certain  applications  it  is  desirable  to  have  analytical 
expressions  for  the  variances.  This  can  be  achieved  by  approximate 
solutions  of  Eqs.  57  fo  60  for  a^  « 1 and  ao  » 1 and  combining 
the  r<  suit  s „ Again,  the  dei'ivat  ions  are  found  in  App.  C;  for  the 
important  practical  case  n = 2 it  is 

for  3o  « 1 

=■■  y ^ [ ( 3 . 8 5 - 3 ^ ao  ) cos®  ePo  + ( 1 . 3 1 - ao  ) sin®  cpQ  ] 
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[Eq,  64a] 


V STl  ^ sin^Cfb) 


_ 8(1  - 2Wo  ) ct  coa>% 


3tt  VP 


(1  -0.596 


= 7.413  • 10"^  cos  cpo  (1  - 2Wo)  (1  - 1.732  y^)  , 

V® 

for  ac  » 1 


°ez  ^ 16(1-  2Wq)  a 

V _ 4 3 

3^0  ^ V ^ 


1 + cos  CPo 


+ 4 00s  Cfb  (-  1 + ao  1 sin  cfe  | ) 
1^2  ^n(2ao)j  , if  cCb  = 0. 


The  dimension  of  a and  radians,  the  other  quantities 

ai’e  dimensionless. 


With  a simple  intorpolat  ion  t'orniula  it  rs  possible  tij  combine  Ih.e 
limiting  cases  in  Eq . 04  to  general  valid  formulae  that  can  more 
easily  be  used  than  th<'  numerical  integration  procedure  described 
earlier  in  this  chapter.  The  interpolation  formula  appropriate 
for  Eq . 64  is 


€(x) 


[Eq.  65] 


The  functions  fi(x)  and  f2(x)  are  the  approximations  for 
a^  <<  1 and  ao  » 1 Eq,  04^  respectively,  with  one  modifica — 

tion:  the  approximation  for  should  read 

“ J ~ [ ( 3 . 8 5 + I a I ) cos^  cp  + (1.31+  a I ) sin^  CP  ] 

' * = I ol  0 


to  provide  a defined  function  for  all  values  . The  exponent  m 
in  Eq.  6 5 can  be  chosen  arbitrarily;  m •■=  2 has  shown  to  approx- 

imate the  equations  with  sufficient  accuracy. 


The  relations  in  Eq.  64 
variable.  For  practical 
chosen.  The  caiculatioii 


contain  the  facet  length  L 
applications  this  length 


+- 


as  the  independent 
L has  to  be 
ne  xt  ch  a pt  e r , 


3.3  Calculation  of  the  facet  length 

The  facet  length  L is  calculated  following  the  general  procedure  of 
Sect.  2.3.  The  starting  equation  is  Eq,  38,  v;hich  is  evaluated  for 
the  Piei'son-Moskowitz  spectrum  [Eq.  55j.  First,  the  variance 

= E^jp^  I is  calculated,  employing  Eq,  42 


2av* 

px  g^ 


TT 

J'  cos‘'(cp  - ePp  ) 


o 


[Eq.  67] 


1 1 

4 u| 


2cos  u 

a. 


4- 


2 


dep 


3 


where  u is  defined  in  Eq,  6I.  The  dei'ivation  of  Eq.  67  can  be 
found  in^App.  C.  From  Eq.  67  it  is  seen  that 


[Eq.  68] 


therefore  the  results  in  Figs.  13  and  14  are  given  in  this  form. 
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As  foi’  Ihc  facet  slat  istics.  the  limiting  forms  for  a^  « 1 aiui 
» 1 arc  calculated  in  App . C.  It  is,  for  k\  = 2: 

for  a„  <<.  1 
0 


ex  ^ - 2 

( 3 cos 

288 


CP  + 


sin^  CPjj) 


for  » 1 

0 


t'"  B 


2g 


[Kq.  69] 


Foi’  tl\csc  limiting  foi’ins,  Eq.  3^  is  solved  in  closed  form,  giving 


for  a.  « 1 
0 


L 1 

^0  '(  3cos'^  cpQ  + sin^cp^) 

[Eq.  70] 

and  for  a,  » 1 


The  two  parts  of  Eq,  70  are  interpolated  to  a general  formula 
with  the  I’esult 


L 


X 


0 


(3 


1 

cos  + sin^epo) 


”4  ? B 


[Eq.  71] 


This  is  an  important  I’osult,  as  it  relates  the  surface  parameters 
of  wind  speed  v aixd  wind  direction  cPo  and  tlie  acoustic  wave 
length  X to  tl\e  facet  length  L.  Equation  7I  is  plotted  in 
Fig.  15,  together  with  the  exact  result,  obtained  by  the  solution 
of  Eq.  38  with  an  iterative  procedure  on  the  digital  computer. 

The  figure  shows  that  the  maxinuun  error  does  not  exceed  20^j  thus 
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the  appi  oxiiiiat ion  should  be  sufl'icient  in  most  cases.  It  is  furtlici’ 
seen  that  the  quality  unnibei'  N — ‘ -L-  [liq.  3o]  is  alwa>s  yreatei’  than 

10,  thei'ei'oi’e  tile  ereoi’s  induced  by°t  he  facets  ace  negligible  and 
t lie  problem  of  choosing  the  facet  length  does  not  need  furthci’ 
im'Cst  igat  ion. 

The  same  procedure  ivas  pen-fornu'd  for  the  exponent  n = 4 of  the 
directivity  law,  the  I'osult  is  shown  in  I’ig.  10.  Figures  15  and  10 
arc  similar,  the  only  difference  is  a larger  variai  i<m  of  L with 
t lie  direction  angle  cpg  . 


3 . 4 Resulting  dopplei'  spect  ra 

Tile  compulei'  results  for  dopplcr  spectra  v^'ilJ  now  be  discussed,  the 
facet  statistics  being  calculated  as  described  in  the  preceding 
section.  The  dopplei'-densit y spectra  use  Eqs.  47  and  48;  the  sea 
surface  spectrum  cp)  emploj'cd  ohcre  is  the  spectrum  due  to 

Scott  [ Eq . A.oO]  and  the  mixing  function  W(f,  tp)  is  found  in  Eq . A.o4. 
This  choice  of  t lie  sea-suri’ace  .spectra  seems  to  be  the  best  possibil- 
ity at  present.  It  should  be  understood,  however,  that  the  general 
t'esult  foT’  Ihe  dopple:  density  of  backseat  tei-ed  sound  can  be  eval- 

uated foi  aii3’  sea-sui'face  spectrum  wi’itta'n  in  a form  used  here. 

In  Fig.  J7  the  influence  of  the  grazing  angle  Yo  the  dopplei’ 

density'  is  demonst  rated.  The  x-axis  shows  the  dopplcr  shift  f - fg 
in  hei’tz,  the  3-axis  shows  the  dopplcr  densit>  cp(  f ) in  l/Hz,  The 
angle  Yo  changed  from  2°  to  20'^,  the  wind  speed  is  v — 8 m/s, 

the  incident  acoustic  frequenej’  is  3.5  kHz,  and  the  main  direction 
of  the  surface  waves  points  awaj'  from  the  sound  source.  This  causes 
an  asymmetry'  in  t lie  spectrum,  which  is  explained  as  follows:  Tlie 

correlation  between  t lie  facet  slope  and  the  facet  vertical  velocity 
is  negat  i\e  in  this  case;  this  means  t Hat  a lai'ge  angle  belongs  to 
a negative  velocitj'.  This  veloci  t.v  cruses  negatiie  dopplcr  shift; 
as  the  slope  is  greater  for  the  negative  velocity  than  foi’  the 
pusitive  velocity’  the  backseat  tored  enexgy  is  great  er'.  In  this 
figure  the  sea-surface  spectrum  is  oiiini  di  lect  j onal  in  the  freqiiencj' 
I’ange  where  resonant  scattering  occurs:  about  2,7  Hz  or  20  cm  wave-- 

leiigl  li.  The  lack  of  s.vrame t r>-  is  due  01  ly  to  the  facet  movement  . 

J t the  orientation  of  the  sea  surface  specturm  is  90®  off  the 
incident  wa\e  vcctoi’,  this  effect  does  not  occur,  as  is  seen  fi'om 
Fig.  IS. 

At  this  point  it  is  instructive  to  compare  these  I’esults  witli  those 
of  another,  i eri'  simple  approximation.  In  tliis  approximation  the 
vai’iances  of  and  arc  calculated  as  second-order  moments 

of  the  sea-surface  spectrum.  This  is  the  special  case  of  the  facet 
model  with  facet  length  L = 0.  The  slope  of  the.  sea  surface  is  not 
considered. 

Figure  19  shows  the  lesult  of  ih.e  approximation  for  a limiting 
frequency  fg  = 1 Hz,  Compared  to  the  previous  figure  it  can  be 
said  t hat  for  small  grazing  angles,  Yo  ^ differences 

are  up  to  20  dB . Mor  eover,  t he  asymmetry  of  the  curves,  wiiich  must 
be  expected  from  phi’sical  reasons,  is  not  incoi'porat  ed  in  the  model. 
The  small  asynm>etries  in  the  fi’equency  region  from  0 to  2 Hz  are  due 
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bocai'o'i  ihe  Bra^g  frequency  is  shifted  to  tlie  onmidirecl ional  part 
of  the  sea-surf  ace  spectrunu  The  l-Liiiitiriy  frequency  is  fg  = S 
an  this  case;  this  means  tliat  for  fo  = 33  Idiz  the  Bragg  frequency 
is  £o  — 8.25  Hz  am)  accordingly  the  spectral  density  cp(f)  is 
almost  syinjiietrical , In  this  fagur-e  the  modified  Scott  sea-surface 
spoctraim  as  used.  If,  instead,  t.lie  unmodified  spectrum  is  used  up 
to  fn  = 5 Hz  and  above  a onanidirect ional  spectrum,  the  i^esult  of 
Fig.  28  is  obtaiiuid.  T)ie  asynunctry  of  the  lower  acoustic  fi'cquencies 
is  considerably  higher  tlian  in  the  previous  figure.  For  fj,  = 33  IcHz 
i lu:  cui'v'os  ai'c  almost  identi  cal.  This  shows  that  the  clioice  of  the 
limxting  fi'cqucncy  and  the  modif icaiioti  of  the  Scott  spectrum  is 
impcrtajii  . A compai'ison  with  measurements  will  be  necessary  to  find 
out  the  optimum  values  foi'  tdiis  jaarameter. 

The  last  tliree  figaiies  show  t lie  I’esults  foi'  the  doppler  density  at 
higher  wind  speeds.  The  limiting  frequency  was  chosen  as  fg=  2,5'Hz  . 
In  Fig,  29  all  pai'amei.ers  arc  the  same  as  in  Fig,  28,  except  that  the 
wind  speed  is  doubled.  In  Fig,  30  tlie  direction  angle  is  changed  to 
ePo  = 90^,  thus  giving  symmetric  curves  in  this  crosswind  case, 

I'igure  3I-  sliows  an  example  for  an  extreme  wi»id  speed  v = 32  m/s, 
where  cA'cn  the  low  acoustic  frequency  fg  = 100  Hz  pnoduces  a doppler 
spectrum  with  a double  gaussian  shape  instead  of  6— function  peaks. 


CONCLUSION 


A 1 hcorc*tical  mode],  foi’  the  doppler  spread  of 
acoUiSi  ic  waves  1 roni  iho  rougli  sea  sm’face  is 
of  a two-compr)uent  str-ucture  (facet  model)  of 


tlie  backseat  tered 
developed  ti'  makes 
the  sea  sun  ace. 


use 


The  model  requires  a description  of  tlie  statisf  ical  pi’operties  of  the 
sea-surface  ioughn>.^ss  including  the  directivity  and  Idle  mixing  of 
incoming  and  out  going  w'avos  i n one  direct  ion  as  a fund  ion  of  f reqiiency^ 
Tlie  concept  i.>f  the- mixing  function  is  necessary  to  obtain  asymmet- 
rical doppler  spectra  as  a function  of  tlie  W'ind  direction.  The 
parameters  of  this  function  need  experimental  verification;  in  par- 
ticular, the  optimal  value  for  the  limiting  frequency  is  not  yet  known. 


Tht?  gcnei’al  result  of  tills  paper  is  a set  of  foi’mulae  for  both-  tlie 
facet  statistics  and  the  dopplci'  spectrum  in  terms  of  the  sea— surface 
covariance  funct  Ion  or  tlic  sea— surface  specti^um.  The  results  can  be 
used  in  two  ways: 


1,  Mcasui'ed  sea-sui'facc  data  are  used  to  fit  the  pai’ameters 
of  analytical  spectra  such  as  the  Piorson-Moskowit z spectrum  01  the 
Scott  spectrum. 


2.  Measured  spectra  are  inserted  into  tdie  general  formulae. 


The  procedure,  developed  foi’  the  choice  of 
that  the  qualit3  number  N — L/^-o  y wliere 
length,  is  nevei’  smaller  than  10,  thus  the 
finite  facet  length  can  be  neglected. 


the  facet  length  L shows 
Xq  is  the  acoustic  wave- 
apei’iur  c ei’i’or'  due  to 
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only  to  the  different  rcsuHs  of  the  small-scale  roughness  back- 
scattering,  They  depend  linearily  on  the  sea-surface  spectrum  and 
have  nothing  to  do  with  the  fac:et  movements  and  the  correlation 
between  slope  and  vertical  velocity.  This  effect  can  be  quite 
important  if  the  sea-surface  spectrum  is  assumed  to  be  directional 
up  to  high  frequencies, 

Fi.gui  e 20  sliows  the  same  case  as  befoi’e  but  for  fg  -=  5 Hz,  where 
the  asymmetrical  effect  is  clearly  seen.  If  the  facet  slope  is 
cor\si  dered  for  the  simple  approximation,  a considerable  change 
occui  s in  the  icsults.  The  st  andai  <1  deviat  ioi\  of  the  facet 

slope  tends  to  iTifiiiity  when  L appr*'aehes  zero.  To  avoLd  numei'— 
ical  instability  a facet  length  of  t =•  0,0001  m was  taken.  In 
Fig.  21  the  r-esult  is  shown.  As  is  to  be  expected,  the  influence 
of  the  glazing  angle  Yo  smaller  and  the  total  hackseatt  ei’ed 

power  is  mmdi  higher. 

These  considerable  differences  show  the  not'd  for  the  taeet  model  in 
the  doppler -spi'ead  calculation. 

The  I'esult  of  the  facet,  theory  lies  between  the  t'Xti'cmes  of  the  two 
pievious  figiii'es  and  is  shown  in  Fig.  22.  Heie  the  asymmetry  of 
tlic  curves  depends  on  both  the  coi’relat  ion  ol'  t h<'  facet  slope  and 
veloriti-  and  t lu-  asymmetric'  r«'su  1 1 s of  the'  resouatit  ba  ckscat  1 cr  ing . 

I'lie  curv'cs  dc^pemd  quit  e st  v tmgly  tui  t h«'  limiting  ft  equency  fg  , 
which  is  specially  apparent  jl  fg  tends  to  infinity,  leading  to  a 
spectiiun  that  has  a directivity  up  to  t he^  liighest  t fequc'ncles . This 
case  IS  shown  in  F'ig,  23  for  the  same'  parameters  as  in  the  previous 
figui  es.  In  fact,  the  figure  cont  ain.s  only  one  bi  anch  of  the  usual 
t wy^  as  only  one  tt-function  occurs  in  the  Bragg-scat t ei  i ng  case, 

TVii  s is  C/t  couT'sc  a somewhat  academic  example. 

Fig'ui'c;  24  dhows  the  variation  of  the  dopplei  density  with  the 

ot  ientaiii  n angle  cp^  for  a grazing  angle  - 6=  , a wind  speed 

V ■-  8 Ill/s,  and  an  incident,  frequency  f^  =■  3.3  kHz.  The  angle  cf> 

IS  vai  it'd  from  0^  t ci  180°  in  steps  of  20*^,  Of  course,  foi  Cp  = 

and  cp,  1 8O0  the  ant  isyinmet  ri  c result  oc-cairs.  The  influence  of 

cp  leaciu-s  f rom  about  5 dB  variat  ion  at  low  dopplei  frequencies  to 
0 

10  dB  at  about  10  Hz  dopplei'  sliitt  . Again,  tl.e  sea-sui  tac  e spectium 
IS  cons  j tU.'i  ed  ointii  g.i  lect  ioiial  at  tlie  Bragg  irequ<'nc>  of  about  2.7  Hz, 

Figui'es  25  and  26  show  the  influence  of  the  wind  speed  foi'  a fully 
developed  si'a.  The  grazing  angle  y is  , the  or  lentation  angle 
cp  is  0°,  and  the  acoustic  frequency  f is  3.3  kHz.  In  Fig.  25 
the  wind  speed  vai’ies  from  2 m/s  to  32  ni^s  in  geonu'Ci'i cal  progres- 
sion, in  F'ig.  26  from  2 to  20  ni/s  liu«;aiily.  Fow  low  wind  speeds 
the  facet  model  Is  seen  to  be  superfluous,  as  the  dopplei'  density 
consists  \ ir(  rally  of  two  6-func'ti  ons  due  to  the  Rragg  resonant 
scat teiing,  lor  high  wind  speeds  the  backseat tcred  energy  is  con- 
siderable liigher,  but  it  is  spread  over  a higher  frequency  range; 
the  bai-kscattered  energy  at  the  Bragg  frequency  does  not  change 
consi  dei'abl  y . 

The  next  figui  es  show  t iie  dependence  of  tlie  i csults  on  t he  incident 
acoustic  fi equency  f ^ . Figure  27  shows  the  dopplei  density  versus 

the  nor  mall  zed  fi’i^queucy  for  f^  - 0.33,  .,  3.3,  10  and 

33  kHz.  F'oi'  the  low  tiequency  t^  — 33C  Hz  tlie  Rayleigh  condition 
is  fulllilled;  thus  the  dopplcr  spectral  density  consists  of  two 
sharp  peaks.  With  increasing  fi equency  the  iacets  become  more  impor- 
tant and  at  the  same  time  the  asymmetry  of  the  curves  disappears. 
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The  general  results  for  the  statistics  of  the  facets  are  evaluated 
for  a Pierson-Moskowitz  spectrum.  ah  appi’oximate  closed-form 
solution  is  obtained. 

For  the  numerical  evaluation  of  the  general  result  for  the  doppler 
spectral  density,  a more  general  sea— surface  spectrum  due  to  Scott 
is  employed,  wliich  contains  the  Pierson-Moskowitz  spectrum  as  a 
special  case. 

Examples  for  different  wind  speeds,  wind  directions,  grazing  angles, 
and  acoustic  frequencies  are  calculated  via  a computer  program. 
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FACET  SLOPE 

Paramet«r;  Wind  airoction  {p 

Fac*t  length:  L [>11]  ° 

y/inft  cnJk«/4<  v/Tm/cn 
'r---.  • L-  - -J 

Directivity  exponent:  n =2 


Correlation  Slope -Vertical  Velocity 


10“^  lO"'*  10"^  10"^  10’’ 


10 


FIG.  7 CORRELATION  COEFFICIENT  P vs 
Pofam*f#f»;  »•»  Fi®.  3a 
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locity  Sjl 


Correlation  Slope -Verti 
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Roughness  on  the  Facet 


FIG.  15  RELATIVE  FACET  LENGTH 

Parameter;  Wind  direction  C(?  3 1 

j Exact  result 
; Approximation 

Facet  length:  L (]m  3 2 

Wind  speed;  v[m/s3 
Acoustic  wavelength;  m] 

Directivity  exponent:  n = 3 

RELATIVE  FACET  LENGTH 

Parameter:  Wind  direction 
Facet  length:  L [^m  j 
Wind  speed;  v [|m/sl 
Acoustic  wavelength:  C'””l 

Directivity  exponent;  n -■  * 


1:  Exact  result 
2:  Approximation 
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FIG.  25 
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f-fo 

DOPPLER  SPECTRAL  DENSITY  ep  (f)  v*  DOPPLER  SHIFT  f - 

Parameter:  Wind  speed;  v«2  m s to  32  m.'^s 
Centre  frequency:  f^”  3.5  kHr 
Grazing  angle:  i° 

Wind  direefion:  (downwind) 

Limiting  frequency:  („  » I Hz 


Kz 


Mo 

fo 


DOPPLER  SPECTRAL  DENSITY  «p  (f)  v»  NORMALIZED  DOPPLER  SHIFT 
Parameter;  Centre  frequency:  fo*  0.3  kHi  to  33  kHi 
Grazing  Ckigle; 

Wind  »peed:  y«#  m/ » 

Wind  direction;  (p<*0“  (downwind) 

Limiting  frequency:  fg»3  Hz 


fo=0-3  kHz 


.002 


FIG.  28  DOPPLER  SPECTRAL  DENSITY  <P  (f)  vi  NORMALIZED  DOPPLER  SHIFT 
Porometert;  tee  Fig.  37 

except;  Limiting  fre^ency;  fg"S  Hz,  unmodified  Scott  spectrum 
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Doppler  Spectrum  tp(f)  ^ Doppler  Spectrum  tp(f) 


-0.008  -0.004  0 .004  .008 


F!G.  31 


I 

DOPPLER  SPECTRAL  DENSITY  ^(f)  v*  NORMALIZED  DOPPLER  SHIFT  - ~ ^ 

Paramet«f$:  $e«  Fifl-  27 

except:  Wind  speed:  v ••  32  m/s 

Limiting  frequency:  fg  «.  2.J  Hz 
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APPENDIX  A 


SEA-SURFACE  DESCRIPTION 


A . I Intro  duct-pry 

The  sea  surface  h(x,  y,  t)  is  considered  as  a stationary  three- 
dimensional  process  that  is  assumed  to  be  gaussian  distributed. 
Then  the  statistical  properties  of  this  process  are  completely 
described  by  its  covariance  function  C , which  is  defined  as 


^ ^ ) Yi  j f 1 > j Ys  i bg  ) E j h ( Xj  j Yi  > bj  ) h ( X2  j Ya  > bg  )!•  A.l3 

The  symbol  e| ...  | denotes  expectation. 

Because  of  the  assumption  of  stationarity  in  space  and  time,  the 
covai’iance  function  is  dependent  only  on  the  difference  of  the 
variables.  The  following  notation  is  iised 


5^1  - Xg  = X,  yi  - yg  = y,  ti  - tg  = T , [ Eq . A . 2 

where  x and  y are  spatial  displacements,  T is  the  time 
displacement . 

The  covariance  function  is  related  to  the  three-dimensional  wave 
parameter  frequency  spectrum  Xa(kx,  k„,  f)  via  the  three- 
dimensional  Fourier  transfoi'm 


+* 


Xa(k^,k^,f)  ==  II  C(x,y,T)  exp(-j2Tr(-kxX  - kyy  + fT)jdxdydT'\ 


and 


> [Eq.  A. 3] 


+» 


rr 


C(x,y,T)  =-  jj  X3  (k^,  ky,  f ) expl  j2TT(-k^x  - k^y  + fT  ) I dk^dk^^df , y 


where  both,  X3  and  C are  real-valued  functions  and  f is  the 
frequency.  The  components  k arid  k of  the  wave  parameter  k 
are  * 


k^  + k^  = 
X y 


and  k 


1 

X 


[Eq.  A. 4] 


where  X 


is  the  wavelength  of 


the  corresponding  surface  wave. 
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Equation  A, 3 states  that  the  sea  surface,  with  the  assumptions  made, 
is  described  by  the  covariance  function  C or  the  wave  parameter 
frequency  function  Xj , 

In  the  following,  both  forms  are  considered  in  more  detail. 


A.  2 The  wave  pai'ameter  fi’equency  specci-um 

Assuming  that  the  dispersion  relation,  (Eq.  4 of  the  main  text)  is 
valid,  the  three-dimensional  spectrum  Xj(k^,k  ,f)  can  be  written 
as  ^ 

Xs  [6(f-f*)  + 5(f+f*)]  , 

y / y y 

[Eq.  A. 5] 


where 


f* 


[Eq.  A. 6] 


In  Eq.  A, 5,  the  function  Xg (k  ,ky)  is  the  wave  parameter  spectrum. 
Xg(k  ,k  ) contains  the  orientational  information  of  the  fi’ozen  sea 
surface^'  The  "mixing  function"  W describes  the  ratio  of  incoming 
and  outgoing  waves  for  a given  wave  parameter  vector.  The  function 
W(x,y)  has  the  properties 

0 < W(x,y)  £ 1 

and 

W(x,y)  = l_W(-x,-y)  . 

It  is  used  in  Eq.  A,  5 with  the  complicated  argument  W[kxSgn ( f ),kySgn (f)] 

to  fulfill  the  symmetry  relations  imposed  on  the  wave  parameter 

specti’a  Xj(k  ,k  ,f)  and  Xg(k  , k ) , namely 
X y X y 


[Eq.  A. 7] 


Xa(k^,k^^,f) 

and 


X.(-k^,  -ky) 


[Eq.  A. 8] 


which  are  necessary  to  obtain  a real  covariance  function. 
Apparently,  the  noi’mali  zat ion  for  Xj  and  Xg  has  to  be 


C(0,0,0) 

and 

C(0,0,0) 


[Eq.  A.9j 
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me  Oasic  equation  LEq.  A.5J  will  now  be  formulated  in  terms  of  other 
functions  describing  the  statistics  of  the  sea  surface  as  well. 

First,  a description  in  polar  coordinai es  is  given.  Equation  A. 9 is 
rewi'itten  as 


= 


Xg  ( k , k ) dk  dk 


X-  y X V 


+TT 


Xg  ( k COS  cp  , k sin  cp)  kdkd  cp  , 


[Eq.  A. 10] 
[Eq.  A. 11] 


o -n 


where  a is  the  standard  deviation  of  the  sea-surface  height  and 
h 


k^  = k COS  cp  , k^  ==  k sin  cp  , 


[Eq.  A. 12] 


Instead  of  the  wave—  parameter  angle  description,  a frequency— angle 
description  can  be  used.  Defining  a function  Fg(f,  cp)  via  the 
normalizing  equation 

* +TT 

= JJ  Fg  (f,  cp)  dfdcp  [Eq.  A. 13] 

o -n 

the  relationship  to  Xgik^jk^)  is 


Fs  (f,  cp)  = 


8tt®  f®  / 2rrf= 


g 


2rrf“' 

cos  cp  , ' '■  sin  cp 

g 


or 


g^  / * K \ 


where  again  the  dispersion  relation  of  Eq.  A, 6 is  used, 
Fg  (f,  cp  ) has  the  symmetry  property 


[Eq.  A, 14] 


[Eq.  A. 15] 


Fg  (f,  cp  ) = Fg  (f,  cp  t rr) 


[Eq.  A. 16] 


Inserting  Eq 

Xa(h^,ky,f) 


A. 15  into  Eq.  A. 5 gives 

" ai’ckaiTj^  )w[k^sgn(f ),  k^sgn(f)] 


[6(f-f->^)  + 5(f+f*)]. 

« 


[Eq.  A. 17] 
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The  mixing  funct-ion  W(k^,ky)  can  be  replaced  by  another  function 
W(f,  cp)  with  the  arguments'’^  f and  cp  . This  function  has  the 
property 


0 < W(f,  cp)<  1 I 

and  I [F,q.  A.l8] 

W(f,  cp)  = l-  W(f,  cp+rr).  ] 


Although  the  argument  of  W has  changed,  the  letter  "W"  is  used  for 
both  functions,  as  there  is  no  danger  of  confusion.  Introducing 
Eq.  A, 18  into  Eq,  A.l?  yields 

^3  f arctan  ^ arctan  |^  + rm(-f)^ 

[6(f  - f*)  + 6(f  + f^)]  , [Eq.  A.19] 

where  u(x)  is  the  unit  step  function 

{0  X < 0 

5 X = 0 [Eq.  A. 20] 

1 X > 0 . 


In  many  practical  cases  the  dependance  of  Fg(f,  cp)  on  cp  is  not 
known  and  therefore  neglected.  Then  a wave-parameter  spectrum 
Xi (k)  normalized  via 


di 


Xi  (k)  dk 


[Eq.  A. 21] 


describes  the  sea  surface.  The  relation  to  Xn ( k ,k  ) is 

® X V 


Xi  ( y k®  + k®'  ) = 2tt  yu  + k®  Xg  ( k , k ) 

\ X y J ^ X y 2 ^ X'  y' 


[Eq.  A.  22; 


Usjug  again  the  dispersion  relation  of  Eq.  ,\.6,  a frequency  spectrum 
Fi(f)  introduced  that  is  related  to  Xi(k)  via 


where  Fi(f)^  is  normalized  with 

^h  = ,/ (f)  df  . 

o 


[Eq.  A. 23] 


[Eq.  A, 24] 
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As 


+rr 

J Fj  (f,cp)  dcp  = l\  (f ) , 
-r 


[Eq.  A. 25] 

i 


the  I'elation  to  Fg  is 


Fi(f)  = 2TTF3(f,w)  , 


[Eq.  A. 26] 


if  F^  is  not  dependent  on  cp  . 


In  the  case  of  omnidi rect tonality  the  mixing  function  W is 
normally  set  to  0.5.  Inserting  now  Eq.  A.  26  into  Eq,  A. 19.  one 
obtains,  for  this  special  case  of  practical  importance 


Xs  (k^,ky,f) 


4(2nf^f)* 


Fi  ( f^" )[  ^ ( ) + 5 ( t'+f*  ) ] . 


[Eq.  A. 27] 


In  this  special 

Y /"  1-  I,  -f  I 

-■3  V ‘'XJ  ‘'■yj  ^ / 
holds,  which  is 


case  the  relation 


not  valid  in  general. 


[ F,q  . A , 2 8 ] 


A , 3 The  covariance  function 

A. 3.1  General 


* -I  T 1 _ j_  j ,.1  , 4 w _ ^ ^ ^ A T *-i  4-  1 vw-i  1 ? 4,1  __  »»  f*  »»  rti  ri  ■{'  f*  ca  n 1 1 ti  r>  xr 

j_  j_  1 • f..  ji.  di  VAC?  iJ.vc;vA  .i_ii  V--I  xca  jj*  t .» j,,  v*  * ^ w«  — - ^ w ^7 

domain  have  consequences  for  the  covariance  function  C(x,y,T) 
according  to  the  Wiener  tlieorcni  of  Eq,  A.  3.  In  this  chapter  the 
relations  are  listed  with  a few  comments, 

Intr'oducing  Eq.  A. 5 into  Eq.  A.3j  one  obtains 


C(x,y,T)  = (k^,ky)[W(kx,ky)  expl  j2rr(Tf*-kxX-k^^y)  1 

+ ^1  - W(k^,  ky)j  expl  j2Tr(-Tf*-kxX-kyy)  } ] dk^  dky  . 


[Eq.  A. 29] 


If  T - 0 , the  ■ ovariance  function  C(x,y,o)  of  the  "trozen" 
sea  surface  is  obtained,  which,  of  course,  is  independent  of  the 
mixing  function  W . 
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Eqvxation  A.  29  can  be  written  in  real  form  using  the  symmetry 
relations  of  Xg  and  To  show  this,  the  following  theorem 

is  applied.  If  f(x,y)  = fj,(x,y)+j  fj(x,y)  = f*(-x,  -y)  , 

(*  : conjugate  complex)  then  its  Fourier  transform  is 


F(X,Y) 


f(x,y)  exp  j 2rrj  (xX  4- yY)  ) 


dxdy 


ZjJ'  f^(x,y)  cos  [2tt(xX  + yY)]  - f^(x,y)  sin[2Tr(xX  + yY)]  dxdy, 
o -» 


[Eq.  A. 30] 


In  this  case  it  is 

f(k^,k^)  =X2(k^,k^^)[W(k^,ky)expl2TTjTf*  } + (I  -W(k^,ky)exp|-2TTjTf*}], 

[Eq.  A. 31] 

which  equals  — k^,  ”^y)  can  be  shown  simply.  Therefore  one 

obtains  from  Eq.  A, 29 

= j j^Xg  (k^,  k^)  { cos  2rrf’''^T  cos  2 TT(k^x  + k^y) 


+ r 2W<"  U u i_ll  ci«2tT  -P*T  o ■:  r,  •>  TT  I-  IT  4-  1.  .ri  1 H L-  H L- 

, , x“*'y 


■ '■“x"  “y'  ‘ 


[Eq.  A. 32] 


Tlie  introduction  of  polar  coordinates  according  to  Eq.  A, 12  gives, 
in  the  complex  form 

*1.  +TT 

C(x,y,T)  =r  j^Xg  (k  cos  Cp,  k sin  cp)  | W(  k,  Qj)  exp[  2it  j ( t f*-k  x-k  y)] 

J J ^ y 

O -TT 

+ [1  - W(k,  cp)]  exp[2TTj(-Tf*-k  X-k  y)]  ] k dk  dq?  , [ Eq  . A. 33] 

y 

where 

W(k,cp)  < 

and 

W(k,cp)  = 


1 - W(k,  cp  tt)  . 


[Eq.  A. 34] 
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The  real  form  is 


■ +n 

C(x,y,T)  = J Jxs(k  cos  cp,  k sin  cp)  j cos  2tt  t f^cos  2tt  k (xcos  CP  + ysin  cp) 
o -n 


+ [ 2W(  k, cp)  -1]  sin  2ttt  sin  2TTk  ( xcos  cp  H- ysin  cp)  j k dk  dcp. 

[Eq.  A. 35] 


Using  again  the  dispersion  relation,  one  obtains  for  the  complex 
form  in  Eq.  A. 33 


+TT 


r p 

^ f 2TTf^ 

2rrf3  . \ 

J 

X-  cos  cp, 

V g 

sjn  cp  1 

^ / 

0 

1 

/ 

2 

W(  f , cp)  exp  I 2rrj[  if (xcos  cp  + ysin  cp)]  | 


+ (1  - W(f,  cp)  ) expj  2rrj[-Tf  - (xcos  cp  + ysin  cp)]  ! ) -f-"  df  dcp  , 

e / 


[Eq.  A. 36] 


wliere  W(f,  cp)  6as  the  property  of  Eq.  A.  34. 


For  the  real  form  one  obtains 
T ( 2nf"  _ . 2^f" 


C(x,y,T)  ==J  I'  X, 


O -TT  \ 


_ 1 COS  cp  , sxn  cp  XCOS 

vs  g j\  S 


( 2tt  f* ) ^ 

2rrT f COS  — ( xcos  cp  + ysin  cp) 


_L.  r OT»7  /'  4'  .a\  1 *1  nmm  -(^rr  J ^ ^ ^ . 4-  , 


■"''J  g 


[Eq.  A. 37] 


Using  Eq.  A.14»  Eqs.  A. 36  and  A. 37  are  formulated  in  terms  of 
Fs  ( f , Cp)  , yi  elding 


• +TT 

p r 

C(x,y,T)  = I j i- 

o - 


Fg(f,cp)(  W(f  ,cp)expl  2ttj[  Tf  - xcos  cp  + ysin  Cp)  ] 

g 


+ [ 1 - W(f , cp)  ] exp  I 2ttj  [-Tf  - (xcos  cp  + ysin  cp)  ] | j df  dcp  , 


[Eq.  A. 38] 
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and 


TV 

r P / Ojrf 

C(x,y,T)  =J  J Fg  ( f , cp)  j cos  2TTfr  cos  [ (xcos  cp  + ysin  cp)  ] 


O -TT 


+ [[  2W( f , ep)  — 1 ] sin  2rrf T sin  — (xcos  cp  + ysin  cp)]  j dcp  df 

6 


[Eq.  A, 39] 

Because  of  fhe  symmetry  properties  in  Eqs.  A,l6  and  A. 34  the 
integration  over  cp  can  be  performed  from  -tt/2  to  +rr/2  , 
taking  the  result  twice. 


A, 3* 2 Special  case 

Now  a special  case  is  considered,  which  simplifies  the  general 
formulae.  The  simplification  is  applied  to  Eq.  A, 39}  for  the 
other  lelations  quite  similar  formulations  apply.  For  the  mixing 
function  W(f,cp)  the  following  equation  is  assumed 


W(f,cp) 


for  |*4J-cp^|<5 


W 


elsewhere 


[Eq.  A. 40] 


Inserting  in  Eq.  A, 39  yields 

+tt/2 


0(x,y,T)  2 


Fg  ( f , cp)  cos  2iTfT  cos  -^-(xcos  cp  + ysinep)]  dqj 


o -tt/2 

CPw+tt/2 


_j_  / otat  t \ r I e J „ r ( — '^’f  ) _i 


A P 


CP^-tt/2 


J * 

[Eq.  A. 41] 


The  second  term  in  Eq.  A. 41  can  be  written  in  this  form  because  the 
integrand  I has  the  propei’ty  l(cp)  = - l(cp  ±tt)  , 

In  the  extreme  case  Wq=1  and  cp^f=0,  Eq.  A.  41  takes  the  form 


« tt/2 

C(x,y,r)  = 2J  j'  Fg  (f,cp)  cos  [2TTfT  - — ( xcos  cp  + ys  in  cp)  ] dcpdf. 

o -tt/  2 

[Eq.  A. 42] 


This  case  is  often  assumed  in  the  literature,  e.g.  [Ref.  lO], 
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For  an  extremely  aimplified  sea  surface,  namely  a pure  sine  wave 
with  frequency  f , tiavclling  in  the  direction  cp^  , Eq,  A. 42  is 
solved  in  clesad  foim.  In  this- case,  F2(f,Cp)  has.  the  form 


^h 

Fe  = —■  L 6(f  - fy,  cp  - epo)  + 6(f  - fo,  cp  - cp^,  - n)] 

and  C(x,y, t)  is 


[Eq.  A. 43] 


C(x,y, t) 


= (2nfo)^ 

cos[2TTf^jT (xcosep^  + ysin  cp^^)],  rf  ICP^!*^  2 


TT 


cos[2rrt\^T  + 


g 

(2ntyP 

e 


(xcos  cp  „ 4 ysin  Cp^  ) ] .,  if  ^ 


[Eq,  A, 44] 


The  other  extreme  is  - 0.5.  Then  the  resu]t  is 

• 4-rr 

C(x,y,T)  = J'  Fj  ( f , cp)  cos  2rrf  T cos  (xcos  cp ysin  cp)  ] dep  df  , 

O -TT 


[Eq.  A. 45] 


If,  in  addition,  Fg  is  indepeiideiit  oi'  Cp  , Eq,  A, 45  becomes 


C(x,y,T)  = f Fi  (f)  cos  2nfr  [UlLlil  + y®  ] df  , [Eq.  A.4-6] 


o 


where  J (x)  is  tlie  Bessel  function  of  older  zero, 
o 


The  complex  form  of  Eq,  A. 4b  is 


, I , . expj  2rrjTf  I 4- expl-2rr.iT  f } (2rrf)^  . — rr-, 

C(x,y,T)=  Fi(l) [ y x^  4-  y*']  df . 


g 


If  the  two-sided  version  of  Fj^(f)  , 

F(f)  = -2  Fi(|fl)  , 

is  used,  Eq.  A. 47  becomes 

4«* 

C(5f-^y,T)  =-■  r F(f)  exp|2TT.iTfj  J [-“^^^-  4-  ] df  . 

^ g 


LEq.  A. 47] 
[ Eq , A . 4 ^ ] 

[Eq.  A. 49] 
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A. 3. 3 Examples 


For  use  in  tlie  application  of  the  general  thooi’y,  two  versions  of 
sea-surface  spectra  tliat  are  of  practical  importance  are,  described. 
The  first  is  the  well-known  Pierson— Moskowiti-  spectrum  that  has  the 
form 


Fi  ( f ) 

Fa(f>cp)  = cos"(cp-cpJ  , 


n 


o 


[Eq.  A. 50] 


where 


g * 


[Eq.  A. 51] 


and 


= j cos”(cp- Cp^j)  = 2 


pa 


n+1 


(^) 


-TT 


r(n+l) 


[Eq.  A. 52] 


n ^ ^ 

II  • « 


[Eq.  A. 53] 


Equation  A. 53  holus,  if  n is  an  even  number.  The  constants  a and 
B in  Eq,  A. 51  are 


a = 0.0081 

and  [Eq.  A. 54] 

P - 0.74 


In  the  same  equation,  v is  the  wind  speed. 

Foi-  calculation  purposes,  bq.  A, 51  is  too  complicated,  therefore 
it  is  slight] y modified  to  give  a spectrum  very  similar  to  the 
Phillips  spectrum.  The  modification  is 


Fi(f)  - 


ag^ 


[Eq.  A. 55] 


0 , else 

wliere  fg  is  tlie  frequency,  for  which  the  equation 


r cLg- 


J 


(2n)  f5 


! o ^ ® I 

expi-0  (-=— r-)  } df  = 


J (2TT)*f6 


df 


[Eq.  A. 56] 


holds.  The  solution  of  Eq.  A, 56  is 
2rrv 


[Eq.  A. 57] 
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The  second  spectrum,  due  t-o  Scott  [Ref.  ISl,  takes  into  consideration 
a frequency-dependent  direct iv it of  the  surface  waves.  The  spectrum 
has  the  form 


(f,cp)  = (f)|cosP  ^ , [Eq.  A. 58] 

where 

S = [Eq,  A.  59] 

2TTf 

and  A^  is  a normalizing  factoi'. 


Unfortunately,  Eq.  A, 58  violates  Eq. 
covariance  function  when  transloi'mcd 
Eq.  A.  58  is  interpi'cted  as  a product 
and  a mixing  function  W(f,cp): 

[,  cp  - cp 
I cos  ( 


A,lh  and  hence  gives  no  I’eal 
into  the  time  domain.  Therefore 
of  a twt>- dimensional  spectrum 


u'  i'  r 

••  t * 


^ '4'V 


[Eq.  A. 60] 
r Km  . A . 6l  1 


and 


t-TT 

r 
J 

-TT  L 


[Eq.  A. 62] 


. " 'Pq  > I . z'P"  'Po,  v3s 

|COs[—; ^ -j  I + |sin( ^Jt 


d«j 


r(2s  + 1) 

22S+3  + 


Equation  A.6l  is  still  unsatisfactory  as  it  does  not  give  tlie 
omiiidirectivity  of  the  high— f r’eqnency  surface  waves. 


[Eq.  A. 63] 
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Therefore  the  Scott  spectrum  is  modified  by  introducing 


w(f,cp)  = 


2B  + 


B + I cos(- 

”1  / > - 

loos  ^ J 


fp  - epo 


)l 


as 


^)l 


[Eq, 


where 


B 


(•f  P 

S 


[Eq 


and 

fg 

is  a 

limiting  frequency. 

For 

for 

fg 

S 

0 , 

-»  - , 

the  spectrum 
it  gives  the 

is  omnidirectional  for  all  surface 
unmodified  Scott  spectrum. 

A.  64] 


A. 65] 


waves. 
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APPENDIX  B 

SMALL  SCALE  ROUGHNESS  SACKSCATT BRING  AND  EXTENSION  OP  THE  THEORY 


The  1 booty  of  backscatl oiiiig  of  acoustic  waves  from  the  moving  sea 
sntl'ace  is  well  established  for  t-hc  case  of  small-sea  surface 
roughrioJis,  For  sake  of  completeness,  ihe  theory  is  derived  in  tltis 
appendix,  together  with  a possible  extension  of  the  accuracy  of  the 
l esults.  However,  this  extension  is  not  further  investigated  in  the 
main  text. 


B , 1 Existing  theory 

The  geometry  of  Fig.  1 of  the  main  text  is  used. 

Combining  the  lesults  of  Ref,  8,  pp.  T3  and  34  the  pi'essure  of  a 
scattered  plane  wav'c  incident  to  an  apcrt.ure  ab  is  appi^oximat  el  y 
in  the  far  field  ( r -♦  en  ) 

sin  Y Pg 

= j k expj  2ttj  (~k  r + f t ) I 

t”*  0 0 0 


+a/2  +b/2 

f k^[x(a  - cos  Yo  ) ■*■  yS  + 2sin  Y„  li(x,  y,  t ) ] | dx  dy  , 


-a/2  -b/2 


[Eq.  B.l] 


wher'e 


a 

3 

f 

0 

k 

0 

Po  ■ 

and 

h(x,y,t  )- 


cos  Y ’ cos  cp 
cos  Y • sin  cp 

fi'equency  of  incident  plane  wave 

wave  parameter  of  incident  plane  wave  = 

pressure  of  incident  plane  wave 
sea-surface  function  as  used  in  Eq.  A.l. 


Tlie  incident  wave  is  assumed  to  be  in  the  x-z  plane  (cp  = n) 
which  can  be  done  without  loss  of  geneiality.  ^ 
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Inserting  Eq.  B,1  into  Eq.  1 of  the  main  text  yields 


f(C,a,p) 


sin'  Y • P, 


+-  a/2  b/2 

f I II  J ko[xi  (a  - cos  Yq  ) + yi3]i  exp|-2TTj  f^t} 

-00  -a/2  -b/2 

exp|2TTj  kj,[xa  (a  - cos  Yq)  + yaP]}  expj  2nj  f ^ (t+ t ) } 
E ( exp{-4TTj  sin  Y^  • [h(x^  ,y^  ,t)  - h(xg  , yg  , t+  T )]  } ) 


exp{ -2TTfT  j dXi  dXg  dyi  dyg  dr  . 


[Eq.  B.2] 


With  the  substitution.  x^  - Xj,  = x,  y^  — yi  = y and  following 
the  derivation  in  Ref.  8,  p.  59  one  obtains 


|(f,a,p)  = 


sin^Y  • P« 

5 exp{-(4TTk  • sinY.)  C(0, 0,0)1 


+»  a b 


r j r (X-  expf  (4TTk.sin  Y*  C(x.v.t)} 

J,  3 w u 


-•  -a  -b 


expj-ZiTjC  (f  - f ) T - X k^  (a  - cos  Yg  ) ■”  y k^P]  ) dx  dy  dT,  [Eq.  B»3] 


It  is  assumed  that  a»(x(  and  b»fyj  for  all  values  of  x and  y 
where  C(x,y,T)  is  different  from  zero.  Then  the  integration  in 


where  Cfx,y,T;  xs  cixtterent  trom  zero,  then  the  xntegratxon  xn 
Eq.  B,3  can  be  performed  to  infinity  and  jn  and  lyi  are  neglected 

against  1,  It  follows  that  ® ^ 


sin^  Y • pf 

* 0 V 


expl-(4rrkpSxn  Y^,)®  C(0,0,0)  ) 


SSS  expj  (4nk|jSin  Yg)^  C(x,y,T)} 


exp|-2nj[(f  - f^)  T - k^x(a  - cos  Yg  ) - k^yp  ] j dx  dy  dr  . 


[Eq.  B.4] 


For  the  slightly  rough  sea  surface,  it  can  be  assumed  that 

(4TTk^  sin  Yq)®  C(0,0,0)«1 


[Eq.  B.5] 


and  the  exponential  in  Eq.  B,4  is  expanded  in  a power  series  and  only 


58 


th©  first  tw(y  terms  are  considered.  It  is 


expj  (4rrkpSin  Yq)^  C(x,y,T)  } 1 + (4TTk^sin  Ygi  C(x,y,T)  , [Eq.  B.6] 

giving,  for  Eq.  B.4 


sin'  Y P«  k® 

;f,a,p)  = ^ a-  [6(£-f^,k^(a  - cos  y„),  k„p) 


0 " 0 ' 


[Eq.  B.7] 


+ (iTTk^sin  Yo)^  X3  (k„(a  ~ cos  Y^,  k^p,  f-fo)]. 


Equation  B.7  makes  use  of  Eq.  A. 3 and  the  integral  representation 
of  the  6— function.  The  first  term  in  Eq.  B.7  denotes  the  specular 

reflection,  the  second  the  scattering  term.  In  this  case  of  back- 
scattering,  Eq.  B.7  takes  a special  form  by  putting  cp  = tt  and 
Y = Yjj  • Using  the  definition  for  the  backscattering  doppler  density 


cp(f) 


t(f,g,p)r^ 

Pc^b 


[Eq.  B.8] 


the  result  is 


cp(f)  = (4n)^  sin4Ypk4  X3  (-2k^cos  Yq,0,  f-f^^ ) . 


[Eq.  B.9] 


This  is  Eq.  3 of  the  main  text.  When  integrated  over  f , Eq.  B.9 
becomes  the  backscattering  strength  q in  Eq.  10  of  Ref,  5. 


If  the  dispersion  relation  [Eq,  4]  is  considered  valid,  the  spectrum 
Xj  (kjc,  ky,  f ) is  replaced  by  the  two-dimensional  frequency  angle 


uiii 


1?  f -e  — \ 

* 3 f'¥/  • 


Applying  Eq.  A,  19  oiie  obtains  for  Eq.  B.9 


2g^  sin^ 
cp(f)=-  ° 


0 


Fa(f^,0) 


[W(f*,TT)  6(f-fo-fJ) +W(fJ,0)  S(f-f^+fJ)] 

[Eq.  B.IO] 


which  is  Eq.  5 of  the  main  text. 

Equation  B.IO  can  be  simplified  by  the  assumption  of  omnidirectionality 
of  the  wave  parameter  frequency  spectrum  for  the  slightly  rough  sea 
surface.  In  this  case,  W = 0,5  (see  App.  A)  and  the  sea  surface  is 
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dcsci'ibccl  by  the  frequency  specti’um  Fi(f).  Using  Eq.  A, 2 
one  obtains 


cp(f)  = 


4 

sin  Y, 


'/s' 


gnk^ 


3 

COS  Y, 


I)  [6(f-f  -f*) 


This  equation,  wlicu  integrated  over  f gives  the  result  known  from 
Eq.  18  of  Ref.  5. 


B , 2 Extension  of  the  thcoi^y  for  the  doppler  spectrum  of  the 
sma]  1— scale  roughness  sea  sui'face 

The  present  theory,  as  is  seen  from  Eq.  B.IO  gives  two  5— functions 
for  the  backscattered  doppler  spectrum  of  an  incident  monochromatic 
plane  waA'C, 

A possiblb  extension  would  be  to  consider  one  more  term  in  the 
series  expansion  of  Eq.  B.6.  This  third  term  has  the  form 

Tg  = ^(4TTk^  sin  Yq)*  C"^(x,y,T).  [Eq.  B.12] 

Here  only  this  term  is  considet'ed. 

Insci’ting  Eq . B,12  into  Eq.  B,4  gives 

('-it  sin^Y  pl  at)  k^ 

4^^^f,a,3)  = 2 1 exp{-(4TTk^  sin  Y^  ) C(0,0,0) 

H-m 

-^(4nkjj  sin  Y53) 


cxp{-2rrj[  (f-fg  )t  - k^x(a  - cos  Y^,  ) - B ] 1 dx  dy  dT . 

[Eq.  B.13] 


Equation  B„13  has  to  be  added  to  the  first  order  result  Eq . B.7 
to  obtain  a second— ordei'  formula. 


Foi’  the  covariance 
is  assunu!d.  Then 


function, 
C(x, y,T ) 


oimii— direct i onality  and  W — 0.5 
is  I’epresented  by  Eq . A,  49,  which 
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Is  iuse  ted  into  Eq.  B.13,  yielding 


5 

H-  0«  M 

= gJJJJJf(Ci)  F(fa)  expiZTTjfiT}  expiZTTjfsT}  ' 

- “ O ; 

7x^4-  y^'  J /x^~rp  j 

expj-2rrj[(f-f|j)  T - k^x(a-  cos  Yq)  - kj^y  B]  1 dfg  dx  dy  dx, 

[Eq.  B.14] 


where 


G = 


^0  Po 


^ ab 


2r' 


(4TTk^  sin  Yg  ) expj  -(4nk^  sinY^)^  C(0,0,0)}. 


[Eq.  B.15] 


Equation  B.14  can  be  integrattid  over  t in  closed  form.  The  result 
is 


4-« 


= G 


r 


F(fi)F(f2) 


[ (2TTfi  )‘ 

g 


*/x®+y‘ 


(2TTfg 


g 


exp  I 2ttj  k^[x(a-cos  Y^)  + P]  1 6 ( f i +fs -f+f  j,  ) dfi  dfj  dx  dy. 

[Eq.  B.16] 


Next  the  integration  ovei-  fj  or  fg  is  performed.  The  result, 
after  integration  ovei-  fg  , is 


5^^^f,a,p)  = G 


GL 


[[2n(f-f,- 

g 


exp  I 2TTJ  k^[  x(a  - cos  Y^  ) "t  Y P ] ! dfj  dx  dy. 


\ 


[Eq.  B.17] 


Changing  ic  polar  coordinates,  where 

x = R cos  Cp  , 
y — R sin  Cp 
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i 


gLves 


+fl® 


= «/Ji 


hr 

0 


rii=£diK].p-i£4^R] 


-»  O 


»TT 


expl  2nj  R V (ct  - cos  )®+  P®  cos(cp  - t(f  ) }R  dR  dcp  dfj  , 

[Eq.  B.18] 


where 

cos 


tt  - cos  Yq 
v^(a  - cos  Yg  y + 8^' 


and  sin  il?  = 


P 

*/(a  - cos  Yft)®  + V" 


Integration  over  cp  gives  a Bessel  function  again 

^0  n r(2TTfx)^  -I  rc2TT(f-fe-fx)r  l 

5^-^^f,a,0)  = Gj  J F(fx)F(f-f^-fx)Jo[ ^ ^J'^oL S 


— « o 


2tt  JqC  2tt  R J~(a  - cos  Y{,  )^  + ] R cJR  df^  . 


[Eq.  B.19] 


Using  a i-esult  from  Ref,  12,  p.  412  and  Ref.  13,  p.  334 


I = J Jg(at)  J^(bt)  Jg(ct)  t dt 
o 


r 1 - a® 

^ if  I cos  a|  = i 1<  1 


[Eq.  B.20] 


TTbc  sin  A 


0 


2bc 


otherwise 


and  inserting  into  Eq.  B.20  gives  the  final  result,  which  is 
written  down  for  the  backscattei'ing  case  using  the  normalized 
backscaotering  doppler  density  as  defined  in  Eq,  2 


Cp 


==  4k  g^  sin  y exp{-(4Trk  sin  Y- C(0,0,0)j 


H# 


J F(fi)  F(f-.f^-fx) 


df, 


ff(f-fo-fi)^  sinA(fi) 


[Eq.  B.21] 
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The  integration  over  is  performed  only  for  those  values  of 

fi  , where  [ sin  A(fi  ) ] < 1. 


For  cos  A > one  obtains 


cos  A — 


fl  (f-r,-fi)" 

■■  ■ ■ ■ - ■ + - - — ' — 

2(f-.f^-fi)^  2ff 


^ ^0  Ve 
2(TTcfi  r 


[Eq.  B.22] 


The  evaluation  of  Eq.  D.21  seems  possible  only  by  numerical  methods. 
This  is  not  investigated  here. 


63 


APPENDIX  C 


STATISTICAL  DESCRIPTION  OF  THE  FACETS 


C.l  Introductorv 


The  starting  equations  for  tlve  derivations  of  tlie  statistical 
pi’opei'tics  of  the  facets  are  Eqs.  12  to  3 8 of  the  main  texe. 
These  equations  lead  to  the  facet  inclination  Ocf  tlie  vertical 


and  hoiizontal  velocities 


and 


and  the  crossnioments . 


Fui'cher  on  in  this  appendix  the  sea-surface  roughness  on  the 
facet  a|^(L)  rs  calculated.  For  aj^(L)  , Eqs.  39  and  40  are 
the  st  rting  point. 


This  appendix  has  two  sections.  Firstly,  a general  result  in 
t er'iDS  of  the  sea-surface  covai’iance  function  C(x,y,T)  is 
derived  which  is  evaluated  in  the  second  section  for  a Pierson 
Moskowitz  sea-sui’face  spectrum. 


C . 2 General  results 

C.2,1  Formulation  in  terms  of  the  covai'iance  function  C(x,y,T) 
The  facet  is  represented  by  a strai.glit  line 

f (x ,, t ) --  c ( t ) •+  b ( t ) , [Eq.  C.l] 

such  that  the  mean-square  eri’or  e(t)  between  the  surface 
lunction  h(x,  y,  t ) | and  the  facet  equation  f(x,t) 

1/2 

e(t.)  = J Lh(x,y,t)  - e(t)  X - b(t)]^  dx  [Eq.  C.2] 

-L/2 


is  a minimum.  Taking  the  dei’ivative  of  e(t)  with  respect  to 
c(t)  and  b(t)  leads  to  the  equations 


c ( t ) - 


1../  2 

r 

1 » j 

-i/2 


X ll(x,y,t)  dx 


and 

L/2 

b(t  ) = i lr(x,y,t)  dx  . 
1 o 
-1/2 


[Eq.  C.3] 


[Eq.  C.4] 
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The 

here 

text 


veiecities  in  vertical 
for  completeness.  It 
that 


L/2 


and  horizontal  directions  are  repeated 
follows  from  Eqs.  15  and  18  of  the  main 


u^(t) 


1 P 3h(x,y,t) 

1 J 

-L/2 


[Eq.  C.  5] 


and 


t 

Uj(.(t)  = ~ f J } - h(-  ^.,y,ti)]  dti+Ux(tg)  . [Eq.  C.6] 

to 


The  roughness  on  the  facet  p^(t)  is  expressed  by  the  mean- 
square  error  e(t),  normalized  by  the  facet  length  L.  It  is 


p"(t)  = f e(t)  , 

i_i 


which  is  identical  to  Eq.  40  of  the  main  text. 

All  these  quantities  are  stochastic  variables  with  the  mean  values 
zero  and  variances  that  depend  on  the  sea— surface  properties  and 
the  facet  length  L . 

The  calculation  of  the  variances  and  crossvar  3. 13  O O £ 3.  S'  S ST*  ^ L'  vx 

in  detail  for  one  example.  The  facet  slope  €(t)  has  the  variance 


,Vi  ,ti  ) Xg  h{xa  ,ya  ,tg  ) dxj,  dxg  } , 

[Eq.  C.7] 

where  t^  = tg  and  y^  = yg  . 

Taking  the  expectation  E inside  the  integral  and  using  the 
rect— function 


E^ie 


• "&I 


L/2  L/2 

/ ^ 
-L72  -l/2 


f I^i^2 

rect  X = ■< 

I 0,  else 


[Eq.  C.8] 


one  obtains 


]_4.A  PP 

JJ  XiXg  rect  — ^*5ct  -£-  E{h(xi  ,yi  ,ti  ) h(xg  ,yg  ,tg  ) } ilxi  dxg  . 

[Eq.  C,  9 ] 
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From  Eq,  A.l  it  is  seen  that  the  expectation  in  Eq.  C,9  is  the 
covariance  function  C(xj  , , yg  jCj  , tg  ) of  the  sea  surface. 

As  stationarity  in  space  and  time  is  assumed,  C depends  only 
on  the  space  and  time  differences.  Therefore  the  substitutions 
X = Xj  - Xg  , y * Fi  - ya  ~ ii  - fa  made,  which  trans- 

form Eq.  C.9  to 


X+Xg 


^t^®^  ^ L^JJ  ^ rect~ 


rect — C(x,  y,  T ) dx  dx^  , 
L 


[Eq.  C.IPJ 


The  integration  over  Xg  is  performed  by  considering  two  cases 


x>0  and  x<0.  For  x>0 

one  obtains 

L 

Eje=|  ,mjc(x,y,r) 
0 

9 -1 

t " 1 

L/2-X 

dx 

-L/2 

[Eq. 

C.ll] 

_ 3 

2 J 

L 

= ^ /c.(x,y,T) 

3-+  ^ 

L L’  ■ 

1 dx  , 

[Eq. 

>. — 1 

H 

• 

0 

o 


For  X < 0 the  result  is 

o 

. rt  O / m\ 


Combining  Eqs . C.12  and  C.13  gives 
L 

E^|e=  } = ^ J C(x,y,T)  /l-  3!^  +2  tlx  . [Eq.  C.I4] 


It  was  mentioned,  before  that  this  expectation  has  to  be  taken 
for  tj  = tg , and  y^  = yg , viz:,  for  T = y = 0.  Then 
C( y^ '^ ) I q-_y_Q  is  3n  even  function  in  x . Therefore  Eq.  C.14 

is  rewritten  as 


EtlcM  - 


L 


x,y,T), 


4- 


h=0 


3 


dx  . 


[Eq.  C.IS] 


which  is  Eq. 


19  of  the  main  text. 
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The  calculation  of  E^|u^}  and  E^fu®|  is  performed  analogously 
3^  C(x,y,T) 


For  E.  ju^i  the  fact  is  used  th.at 
"t  ^ 


>'--=1=0 


is  an  even  function  in  x,  while  E ju® } is  calculated  using  the 
relation  that  j Uj^ } = 0 and  therefore  C(x,y,r)  = 0 for  t . 
The  resvilts  are  Eqs,  20  and  21. 

For  the  crossmoments  a similar  derivation  applies.  For  j e u^ } 

one  has 


E/2  L/2 

Xi  h ( Xi  ^yijti)  u(  Xg  , yg  j tg  ) dxj  dXg  • [^Eq.  C.l6] 


Et{^-xl=  TT 

-L/2  -L/2 


Again,  the  substitutions  yi  - = y , - tg  = t and  for  the 

f ii’st  summand  x^  - ^ = x , for  the  second  + ^ = x,  at*e  ma 

It  follows 


o T 

E|c  ^ -j  J(x+t)  I C(x,y,T')  dr*  dx 


^.-L 

1.  T 

X - I'j  J-  C(x,y,T  » ) dT  « 


T“0 

y=0 


ly=0 

dx 


[Eq.  C,17] 


As  C(x,  y,  T ) j is  even  in  x,  the  integral  over  C is  uneven 

in  X , Therefoi’e  both  summands  cancel  out  and  the  r'esult  is 


Et ! e j = 0 . 


[Eq.  C.l8] 


The  same  argument  leads  to 


[Eq.  C.19] 


This  is  Eq.  23  of  the  main  text. 

For  E^  { € u I the  result  is  different  from  zeroi 

T ' 'T  * 
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It  is 


L/2  . 

-jj  E|l.(x.,yi,t.) j- Idx, 

~L/2  -L/2  ” 


t 

j 


dxo 


[Eq. 


12  pr,  . X . 

- JJ  (X  + X,  ) roct  ^ 


X + Xj,  Xa  ~ac(x,y,T) 


rect 


St 


dx  dXa 
T = 0 = 

y=0 

[Eq. 


The  integration  over  Xg  is  performed  for  x > 0 and  x ^ 0. 

-L/2-x 

-1./2 


X > 0 that 

L > ^ 

n 9C(x,y,T) 

x| 

J 9t  |t-0 

* ^2  ~ 2 

o (y=0 

_ 

dx  [Eq, 


12  rac(x,y,T) 

L*  .1  3t 
o 


L x^  \ , 

7 x - — dx  . 


r=n' 


[ Eq . 


For  X < 0 the  result  is 


E^ic  r,  } = 


_ 12  P9C(x,y,T) 

J 3t 

-L 


y-T=0 


— X + — I dx  . 

2 2 j 


[ £q> 


The  function 


f(x)  = 


1 X 

-■  X - — , if  X > 0 

2 2 


I — X + — , if  X < 0 
'■2  2 


IS  uneven  in  x . As 


3t 


[Eq. 


is  also  uneven  in  x , the  result 


C.20] 

C.2l] 

It 

C.22] 
C.2  3] 

C.24] 

C.2S] 
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j2  P 9C(x,y,T) 


wl  icli  is  Eq.  2 2 of  llie  main  text  . 

Next  the  roughness  | p*'  ) is  calculated.  Insert  ing 

Eq.  C.2  into  Eq,  C,6  and  taking  the  expectation  gives 


L/2 

L J 

-L/2 


E![l>(x,y,t)  - e(t)  X - b(t)]^i  dx 


[Eq.  C.27] 


Inserting  Eqs . C.3  and  C.4  I'or  e(t)  and  b(t)  gives,  after  a 
straightforward  calcul at  ion. 


\ * 


L/2  L/2  L/2 

Ejh^  jdx  - J Ejxih(xi,yi,ti)x2h(x2,yg,tj)ldXi  dXj 


-L/2 


-L/2  -l/2 


L/2  l/2 

1 p r 


E { h ( X;^  , y.  , t , ) li(x2  , ys  , tg  ) } dXi  dxj 


-l/2  -l/2 


[Eq.  C.28] 


Now  the  integration  technique  used  for  the  facet  statistics  is 
applied  again,  giving 


1 2 r p ^ ^2  ^2 

} = c(0,0,0)-j^-  Jj  (x  + X2)x2rect — rect  — C( x,  y,  T ) dx  dx^ 


pn  X 4 Xs  X3 

— I reel  rect  C(x,y^T)  dx  dXg 

L u J L L 


[Eq.  C.29] 


The  integration  over  Xj  yields  the  final  result 

L 

A.  1 . X_ 


E^i pM  = C(0,0,0)  - J J C(x,y,T)  ^ ^ 1 dx,  [Eq.  C,30] 

^ -o' 

y=0 


wliich  is  Eq.  4^  of  the  main  text. 
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1 

C.2,2  Forniulatiou  in  terms  of  the  frequency  angle  spectrum 
J'a  (f^cp) 

The  dci'i vations  make  use  of  the  relat ionsliip  between  the  covariance 
function  C(x,y,T)  and  the  frequency  angle  sea-surface  spectrum 
Fg(f,cp),  This  relationship  is  calculated  in  Eq.  A, 39  and  is 
lepeated  here  foi’  convenience 


TT 


C(x,y,T)  =f  f Fg  (f,cp)  j cos  2TTfT  cos[-^^--^—  (x  cos  CP  + y sin  cp)] 
u (J  S 


-TT  O 


+ [2W(f,cp)  - 1]  sin  2rrf  T sin  [— — ( x cos  cp  + y sin  cp)  ] | dep  df  . 


S 


[Eq.  C.31] 


This  equation  is  inserted  into  the  expressions  for  the  facet 
statistics  and  the  roughness  on  the  facet  ojj.  Then  the 
integration  over  x is  performed.  The  calculation  is  given 
for  one  example,  all  other  formulae  follow  in  a similar  manner. 

From  Eq.  C.15 


(l-  f ^-^Jdxdcpdf  . 

[Eq.  C.32] 


L n « 


E J I 


M 


If 


Eg  ( cp)  cos 


O -TT  O 


(2TTf)‘ 


X cos  cp 


As  T=0,  the  mixing  function  W(f,cp)  has  no  influence. 
The  integration  over  x is  performed  giving 


E^jeM 


n « 

24  r P 

L®  J ij 


Es  ( f j cp)  { ( 1 + cos  at)  — — t ( 1 - cos  a L ) 


La*^ 


12 

L*a* 


12 

1="  a® 


sin  a L } df  dcp  , 


[Eq.  C.33] 


where  the  abbreviation 

. ==  V-rrf)^ 


cos  cp 


g 


[Eq.  C.34] 


is  used.  Equation  C.33  can  be  shown  to  be 
TT  » 


P P Es(f,cp) 

/,  aL 

La  cos  — 

- 2 sin  — ) ^ 

1 J a* 

V 2 

2 / 

dfdcp  [Eq.  C.35] 


-TT  O 


which  is  Eq.  24  of  the  main  text 


The  calculation  foi’  the  oilier  moments  anti  the  ciossmoniciit  is  not 
made  here  as  it  is  sti’aight  foi'wai  d and  follows  the  above  derivation 
strictly.  The  resulting  formulae  are  found  in  the  main  text. 


C.2.3  Limiting  cases 


To  gain  a bettei'  understanding  of  the  physical  meaning  of  the 
expressions  dei’ived  in  the  preceding  section  it  is  instructive 
to  considei'  two  limiting  cases.  Tliey  ai’c  connected  with  the  cor- 
relation distance  Xj^  of  the  surface  covai’iance  function  C(x,y,T). 
The  distance  xr  is  defined  as  the  5 dB-width  of  C(x,y,r)  in  the 
x-dircction.  The  two  cases  L » Xj^  and  L « xr  are  considered. 

a)  L » Xr 


If  the  facet  length  L >>  Xr,  then  C(x,y,T)  i:s  0 for  all  values 
of  X in  the  order  of  L , that  means,  2i  can  be  neglected  against 

1 in  the  formulae  for  Ihc  facet  statistics.  This  gives  immediately 
the  limiting  cases 


IT  J C(x,y,T)|y.__^^Q  dx 


o 
L 

2 pa^C(x, y,T) 


o 


T ti 


^ dx  , 
y=T-=0  ^ 


Eju^)  fv ^ [*  f C(0,y,r)  dt  dti  , 

X _ C .1  ' iv—  i—ii 


T « g U 

^09 


I > ■“  ' 


[Eq.  C.36] 
[Eq.  C.37] 


[Eq.  C.381 


and 

^t  ^ ^ ^ 


L 

1_2  p3C(x,y,T ) 

J 9t 


y=T=0 


X dx  , 


[Eq.  C.39] 


For  the  roughness  on  the  facet  j } follows 

L 

E^ip=  } C(0,0,0)  - I f C(x,y,T)|^^^^  dx  . [Cq.  C.40] 

o 


From  the  equations  it  is  seen  that  all  facet,  statistics  [ Eqs , C#36 
to  C.39J  tend  to  zero  for  L -»  » while  j p^  } equals  the  sea- 
surface  roughness  C(0,0,0)  in  the  limiting  case  L *• , 

Again,  the  above  expressions  are  foi’mulatcd  in  terms  of  the  frequency 
angle  spectrum  Fg(f,Cp). 
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w 


For  cumiilcLencss  the  results  arc  listed  below 


TT  • 


,(31  24  p p *‘a  ^ 

I I — J J sin  at  dt  dep  , 


L®  J t a 

-TT  O 

TT  “ 


~ I ‘’•(‘■"'I’)  ‘"P  < 


COS  Cp 


-TT  O 
TT  o 


r,  ! 3 1 p p Fa 


L“  U u (2ttC)= 

-TT  O 
TT  » 


[Eq.  C.41] 


[Eq,  C.42] 


[Eq,  C.43] 


/ Fs  (f,cp)  [2W(f,cp)  - 1]  2TTf  dep  , 


-TT  O 


[Eq.  C.44] 


and 


q <• 


TT  « 


1p^  } JJpa  (f,cp)  df  dep-  J Fa  (f,cp)  df  dep  . [Eq.  C.45] 


-TT  O 


-IT  O 


b)  b « 


In  this  case  it  can  be  assumed  tliat  C(x,y,T)  changes  very  little 
in  the  integration  range  0 to  L . Therefore  C(x,y,T)  is 
expanded  in  a power  series  and  the  first  two  terms  are  considered. 
This  gives  the  limiting  cases 


t ' 


fki  — 


c(x,o,o) 

ax=^ 


' x=0 


[Eq.  C.46] 


E , j u^ 
t ^ z 


— 


a^  c(o,o,t) 

ar^ 


T=0 


T t. 


E^{u^i 


a^C(x,y,T) 

ax^ 


dt  dt 


^ X=\'=T=0  * 


a^C(x,y,T) 


ar  dx  [ X=y=T=0  ' 


[Eq.  C.47] 

[Eq.  C.48] 
[Eq.  C.49] 


and 


E^.  ! 0 . 


[Eq.  C.50] 
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In  tcrmci  of  the  frequency  angle  spectrum  above 

equations  are 


Fa(f,cp)  df  dtp  , 

Fa  (f,Cp)  (2nt)^  df  dep  , 


TT  « 

Etlu^l  Fg(f,cp)  (2rrf  cosep)^  df  dep  , 

-rr  o 


and 


rr  « 


E,  leii  } «rr  Fj  (f,cp)  [2W(f,cp)  - 1]  2TTfadfdcp 


-TT  O 


[Eq.  C.5l] 


[Eq.  C.52] 


[Eq.  C.53] 


[Eq.  C.54] 


C . 3 Application  to  a Pierson-Moskowitz  spectrum 


C.3.1  Solution  for  general  directivity 


Eor  pract'j-v^al  puipOi^^'s  geuei'ai  1. oi'ninjae  i. n r*qs.  24  to  2C5  of 

the  main  text  have  to  be  evaluated  for  a special  frequency  angle 
spectrum  Fg(f,cp)  and  a mixing  function  W(f,cp).  The  functions 
chosen  are  a Pierson-Moskowit /,  spectrum  and  a simple  mixing 
function,  wliicli  is  found  in  App.  A in  Eqs.  A. JO,  A. 51  and  A. 40, 

First,  the  expression  for  | C ( is  derived.  Inserting  Eq.  A, 50 
with  Fq.  A, 55  into  Eq.  C.33  gives 


24  An 


rr  • 


If 


cos 


'"(CP  - 


--2-|  (1+cos  at)  + (1  - cos  at) 

L.S  L>~’  a ■* 


-rr 


12 

a= 


sin  al. } df  dep  . 


[Eq.  C.55] 


In  Eq.  C.55,  the  factor  a is  an  abbreviation  seen  from  Eq , C.34 

and 


ct  g^ 

A — ■ y 

(2tt)^x 

n 

viheve  x is  defined  in  Eq , A.  52. 
n ^ 


[Eq.  C.56] 
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Equation  C.55  is  integrated  ovei'  f in  closed  form,  Tlie  derivation 
is  lengthy  but  sti'aigbtfoi’ward,  .1  is  denoted  that 


I.,  = 


f 5 


[■ 


( i + cos  bf^  ) 


3L 

715- 


+ ( 1 - cos  bf' 


12L 
^ b^  f8 


12L 

b»f6 


sinbf^  df 


I 


Feo. 


where 


( 2tt)^  L cos  cp  aL^ 


[Eq.  C.58] 


With  the  substitution  u = bf^  and  continuous  application  of 
partial  integration,  Eq,  C.57  becomes 


3L 

2fs 


+ cos  u. 

0 


+ sin  u 


0 


_1 

36u^ 


»Q  p COS  u 

72  J ~ 


[Eq,  C.59] 


where 

Uq  = ao  cos  cp  = bfg 


and 


a = 
0 


( 2nf L 


g 


^/0  Eg 

“7 


[Eq.  C.6o] 


As  ir*"t *3 ^1? i o Cp  is 

computer,  the  problem  of  numei’ 
Therefore  an  approximation  for 
expanding  the  sine  and  cosine 
and  taking  the  first  few  terms 


T^or»'P#-»T>moH  -iri  Hiir-i'f-ral 

2_r  w «.  w > ^ ^ ^ j ^ ^ 

ical  stability  occurs  for  small  u , 

I u^  I <0.2  is  calculated  bj'^ 
terms  in  Eq.  0,59  in  a power  series 


This  leads  to 


3L 

2fl 


f"  7 u^ 
c 

288 


0 

320 


cos  u 


du 


u 


if  |li  |<0.2  [Eq.  C.6l] 
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Inserting  the  result  for  in  Eq.  C.55  gives  finally 


n 


n.  X f 18  . 48 

cos  (cp_  cp_J  4 


'^0  ' 1 2 - A 

0 0 


J-  48  \ / 2 48  \ 

° V u*  / ® \ 0 u u»  / 

" « '00 


- ‘"o  Ci(|uJ  ) 


dcp  j if  >0.2  y 


[Eq.  C.62] 


and 


a V* 


Px  L=g^ 
f n *= 


-J  cos"(cp-cp^) 


^-^-u®Ci(|uo|  )|  dcp  , 


if  u <0,2 
0 


This  is  Eq.  57  of  the  main  text. 

For  the  other  facet  statistics  the  calculations  are  made  in  a 
similar  way.  The  main  steps  are  listed  below. 

Inspection  of  Eqs,  25  and  26  of  the  main  text  shows  that  E^ju^} 
and  E 5u”  j aic  identical  in  respect  to  the  integration  over  the 

V X 

frequency  f.  Inserting  the  Pierson— Moskowitz  spectrum  into 
Eq.  25  gives 


TT 

2 g^A^  p p cos"(cp  - cp„  ) (1  - cos  at) 


\1<I  - -TT-J 


-TT  f. 


f 5(2TTf  cos  cp)^ 


df  dcp  . [Eq.  C,6  3] 


The  integral  1^.  is  in  this  case 
pi  - cos  bf^ 


— df 
f ®( 2nf  cos  cp)^ 


[Eq.  C.64] 


with  the  solution 
(2ttL)^ 


f 6g^f|u| 


1 + 


. V?  \ 

( -1  + — I cos  u + — sin  u 

V 2 / 0 2 


+ 


_r  0 ' j Sin  u 

0 ^ 2 


du 


u 


\^o\ 


[Eq.  C.65] 
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The  approximation 

of  Eq.  C.65 

for  small 

J.  __ 

j 3 

3n|  i.9u^j 

«8 

- + r 

''  6g2  f3» 

U 

8 720 

^ / 

lu 
• 0 

du  V,  if  I I < 0 . 2 


Tnsei'fing  lJiis  J nto  Eq . C.63  gi\  es  the  result 


[Eq.  C,66] 


TT 


E rcas*^(cp-  Cp  ) I -^  + cos  u 

-VPx^u  0 


1 _I 


0 2 


Sin  u I 1 

+ ~-r—  + Si  (|u^I)  dcp  , if  I I >0.2  , 

0 


[Eq.  C.67] 


and 


2ai-° 


TT 


3 Vs  X 


■J  COS^(((l~  (pj  j ^ 


+ 


19ut, 


+ S r ( I U I ) ^ dec 

720  2 ® 


if  juj  ^ 0.2 
For  E, ! u® } the  result  i 

t.  ‘ X ’ 

E |u^  I 

t ‘ X ’ 


2 r n,  ^ ( 
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[Eq.  C.68] 
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Lquaiions  C.6)  and  C.68  are  Eqs.  58  and  59  of  the  main  text, 
r c s pc  c t i V c 1 }' , 


■o,iL>r  iiipniii  ' 


The  ci  o^smoment  1 S | is  derived  in  the  following  way»  With 
the  Pierson-Moskowit z spectrum  ana  the  mixing  function  W(f,cp) 


W(f,cp) 
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[Eq.  C,69] 


Equation  2'/  of  the  main  text  has  the  form 
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[Eq.  C..70] 


As  W(f,cp)  is  independent  of  f in  this  special  example,  the 
integral  over  f has  the  form 
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The  result  of 


the  integration  is 
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[Fq.  C.72] 


where  S(x)  is  tlie  Fresnel  integral  and  the  constants  ki  to  k 
have  the  values  k;^  — - 10/63  ) ~ 8/189  > ka  = - 4/9> 

= - 4/63  , and  kj  = I6/I89. 

Now  Eqs.  C„72  and  C.6o  are  inserted  into  Lq.  C,J0,  giving  the 
desired  result 
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[Eq.  C.73] 


The  factor  before  the  integral  is 
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[Eq.  C.74] 


then  Eq.  C.73  equals  Eq,  60  of  the  main  text.  I'or  small 
|u^l  < 0.2  , Eq.  C.73  is  approximated  by 
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[Eq.  C.75] 


The  roughness  on  the  facet  j } is  derived  in  a similar  form. 
The  starting  equation  is  En . 43  of  the  main  text,  in  which  the 
Pleison-Moskowitz  spiecti’um  is  insci’ted.  It  is,  in  this  case. 
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[ Eq. 


C.76] 


r?ic  rc.-ult  of  the  integration  is 
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[Eq.  C,77] 


Fcr  small  I I '^0,2  , this  is  expanded  in  a power  series  and 

the  first  terms  are  taken.  Inserting  Eq.  C .77  into  Eq,  43  gives 
the  result  for  I"qI  ^0*^  and  |uoj<0.2 
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[Eq.  C.78] 
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This  is  Eq.  67  of  the  main  text. 

The  factor  befor'e  the  integral  is  the  explicit  expression  of 
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[Eq.  C.79J 


C.3.2  Appr'oxLmale  Analy*  ic  solution  for  cosine  square  directivity 

The  expressions  of  Eqs,  57  to  62  of  the  main  text  derived  in  the 
preceding  section  are  evaluated  ’ere  for  the  case  n = 2,  the 
cosine  square  directivity  law,  to  obtain  an  approximate  analytic 
solution  that  does  not  need  a digital  computer  for  numerical 
evaluation.  For  most  practical  cases  these  approximations  are 
cons idtii'ecl  of  sufficient  accuracy. 

In  this  section,  the  expressions  of  Eq.  64  are  dei'ived, 

C.3.2. 1 - ./e^Ic^  r 

a ) a 0 « 1 
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■^0  g.  •-= 


70 


The  starting  point  is  Eq.  Co62  for  cos  cp  < 0.2, 

The  sine  integral  Ci( iu^ | ) is  expanded  in  a series,  giving 
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([  Eq . C 


The  underlined  tei'ms  in  Eq.  C.80  are  neglected,  with  the  result 
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[Eq.  C 


J'h  1 s integral  can  be  solved  in  closed  form  with  the  substitution 
CO  s cp  — X , the  resul  t being 
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[Eq.  C 


Equation  C.82  equals  a®  in  Eq.  64  of  th ^ main  text. 
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0 

If  a » 1,  this  also  has  the  consequence  that  u = a,  cos  cp 
0 0 6 0 

is  much  greater  than  one  in  most  of  the  range  from  0 to  TT, 


However,  if  cp  = '^  , then  u = 0,  ar 


^ ^ ^ , „ ^ 1 4 4- 
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the  integral  in  Eq,  C.62  into  two  parts,  one  for-  u^  « 1 and 
one  for  » 1,  If  » 1,  the  cosine  integral  Ci(|u^|)  in 
Eq.  C.62  is  replaced  by  its  asymptotic  expansion  for  large 
arguments.  Then  it  can  be  shown  that,  neglecting  all  terms 
containing  powers  more  than  l/u^  , Eq.  C.52  takes  the  simple 
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[Eq.  C 


For  the  integration  range  cp  = ~ — € ’ to  cp  = — t e ' , the 
approximation  for  small  u^  in  Eq.  C.8l  has  to  be  taken.  The 
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sum  of  both  gives  the  desired  result.  Tlie  choice  of  c’  contains 
a c<5rtain  amount  of  arbitrariness;  ifor  lack  of  a better  possibility 
X 

e’  = — was  taken. 

Then  it  can  be  shown  that  the  part,  of  the  integral  where  « 1 

can  be  neglected,  therefore  Eq.  C.83  gives  the  complete  result. 

The  integration  of  Eq.  C.83  is  straightforward,  leading  to  the 
t oi'in 
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[Eq,  C.84] 
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The  tan  cp  tei’ins  ds  expanded  for  arguments  near  y , giving 

the  result 
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TT 


sin^ 


which  is 


part  of  Eq. 


64 
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[Eq.  C.8S] 


a ) « 1 

From  Eq.  C.67  it  is  seen  tliat  for-  a^  « I the  vertical  velocity 
has  the  variance 
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[Eq.  C.86] 


which  leads  easily  to 


z ^ a 


This  is  the  result  given  in  Eq,  64  of  the  main  text, 


[Eq.  C.87] 
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Again,  Eq 


b)  » 1 

. C.67  shows  the  approximation  for  a^  » 1 
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[ Eq ♦ C , 8 8 ] 


as  all  sine  and  co.sine  terms  cancel  out  if  the  sine  integral  is 
expanded  into  an  asymptotic  expansion.  This  is  a form  of  the 
integrand  idenlical  to  that  of  Eq.  C.83.  Therefore  the  same 
argument  applies  and  the  result  is 
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which  is  listed  in  Eq.  64 


[Eq.  C.89] 
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The  asyr.'|>totic  expansion  of  the  sine  integral  si  (x)  in 
Eq,  C,6Q  shows  that  also  here  the  sine  and  ct  ine  terms  cancel 
out  leaving  an  integi’al  oi  the  form 
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The  evaluation  of  Eq . C,90  gives 
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[Eq.  C.91] 
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The  sine  and  cosine  terms  cancel  out  i.n  t!iis  case,  too, 
follows  from  Eq . C.68  that 
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which  leads  to  the  result 
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C.3.2.4  E^jeu^l 


a ) a « 1 
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Foi'  small  aiguinents  the  sine  integral  in  Eq.  C.75  tends  to  zero. 
The  term  with  u®  is  neglected  against  u^  , leaving  an  integral 
c)l  tlie  form 
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For  most  practical  cases  is  seems  reasonable  to  assume  that 

c,  p = cp^^.  For  tlixs  special  case  the  above  integral  is  evaluated. 

Application  of  known  integration  formulae  gives 


24«  ( 1-2  Wq  ) cos  y,,  16  yir 


where 


Y ^ r»(  2^)  2"-i 

’'n  “ J <*»  - . [Eq.  C.96] 

O r(n+l) 

Nixmerical  evaluation  of  Eq.  C.95  and  Eq.  C,96  «ivos  Eq.  65  of 
the  main  text. 


b)  » 1 

This  case  leads  to  a similar  form  to  that  for  the  standard 
deviation  Og  of  the  facet  slope*  Starting  from  Bq.  C.73,  all 
terms  containing  sine  and  cosine  functions  cancel  out,  the 
fresnel  integral  can  also  be  neglected  in  this  case.  The 
remaining  equation  is  simply  for  *» 


<Up  . [Eq.  C.97] 
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But,  as  it  is  known  that 
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as  IS  seen  from  Eq,  C.94;»  again  the  technique  of  splitting  the 
integral  is  performed.  The  contribution  for  small  cos  cp  can  be 
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neglected.  Two  cases  have  to  be  considered 
The  first  case  leads  to 
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[Eq,  C.98] 


The  evaluation  gives  approximately 
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The  general  case  Cp^  0 leads  to  the  form 
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[Eq.  C.99] 
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[Eq.  C.lOO] 
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The  result  of  the  integration  is 
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\vhere  the  approximation 
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is  used,  Eq , C.lOl  is  a part  of  Eq.  64  of  the  main  text. 
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Starting  from  the  second  part  of  Eq . C./S  one  obtains 
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which  leads  directly  to  the  desired  result 
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This  is  part  1 of  Eq,  69  of  the  main  text, 


C.lOl] 


C.102] 


C.I03] 


C.104] 


85 


b)  a » 1 
0 

For  a large  > the  bracket  in  Eq.  C.78  (first  part) 

to  as  all  other  terms  can  be  neglected  against  it, 
integral  is  solved  quite  easily.  Starting  from 

S TT 

gives 

X 

^ ^ 2g 

which  is  part  2 of  Eq . 69  of  the  main  text. 
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[Eq.  C.106] 


APPENDIX  D 


CALCULATION  OF  THE  DOPPLER  SPECTRUM  FOR 
AN  ARBITRARY  SEA-SURFACE  ROUGHNESS 


The  starting  point  for  the  calculation  Is  Eq.  46  which  is 
repeated  here  for  convenience 
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[Eq.  D.l] 


Making  use  of  Eq,  6 and  replacing  by  — gives  a form 

of  Eq,  D.l  better  suited  for  the  following  calculations 
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[Eq.  D.2] 


To  solve  this  threefold  integral,  Eqs.  44,  45  and  29  are  inserted 
and  the  integrations  are  yierformed  step  by  step.  As  it  is 
assumed  that  a is  a siiiaJl  angle,  the  following  approximation 
J.S  valid  ® 
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[Eq.  D.3] 
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rhf  ijfst  integration  is  porfoimcd  over  t.lie  variable  u^, 
Because  of  the  5~tunction,  this  is  possible  by  setting 


r - ± = 0 


[Eq.  D.4] 


atid  inserting  irUo  Eq.  D,1  the  value  ^ that  fulfils  Eq.  D.4. 

It  is 
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[Eq.  D.5] 


Inserting  Eq.  4S  into  Eq.  D<,5,  tiie  square  root  can  be  drawn 
appi-o^imately^  as  e and  a , are  very  muclr  smaller  than  c . 

Applying  Eq.  D.2  and  using  the  abbreviation  Eq.  6 
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the  i-esult  is 
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[Eq.  D.7] 


As  f*  is  very  much  smalJei'  tlian  it.  is  neglected  in  the 

tei'iiis  with  and  u^,  Furthcrmoi’o,  sin  yp  replaced  in 

Eq  r D.7  by  siny^  ,«  which  simplitics  the  mathemat  i cat  effort 
considerably  and  causes  a negligible  error  in  all  practical  cases, 
where  is  a small  ang'e.  Then  Eq . D.7  becomes 
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[Eq.  D.8] 


Insert  ing  it  into  Eq.  45  gives,  foi‘  the  doppler  frequency  fp  , 
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[Eq.  D.IO] 


Insei“Cing  ^*<1**^*  D.8  and.  D,9  into  Eq.  D«2  and  using  Eq*  gives 
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[Eq.  D.12] 
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The  function  W2(e,  u ,)  is  rewritten  in  the  form 
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[Eq.  D.I5] 


Using  tlic  well  known  I'olai  ions  for  the  higher  momenfs  oi  i ;i<;  noi'inaJ 
dislribuUion  density,  Eq,  D.13  is  integrated  over  € , yielding 
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[Eq.  D.I6] 


Tlie  product 
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JS  split  up  into  a term  containing  t.lie  variable  ar.d  one 

without  u , , It  is 
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where 


( f - f 4 r ) c 
j 0 0 

2 f sin  Y 

0 ' 0 


[Eq.  D.19] 
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[Eq.  D.21] 
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[Eq.  D.22] 


To  obtain  the  final  result,  the  integration  ovex’  u^  is  now 
pei'formed,  applying  again  the  relation  for  the  moment  calculation 
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[Eq.  D.24J 


and 
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+ tan  y 


[Eq.  D.25] 


The  combination  oC  Eqs.  D.12  and  D.23  g.Wes  the  general  result  of 
Eqs.  47  and  49  of  the  main  text. 
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